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C^ 1 Abstract 
& 

We present a detailed study of AF = 2 observables and of rare K + {Kl) and B s d 
meson decays in a "Minimal Theory of Fermion Masses" (MTFM). In this theory 

00 ■ Yukawa couplings are generated through the mixing with heavy vectorlike (VF) 

fermions. This implies corrections to the SM quark couplings to W^, Z° and 
■ Higgs so that FCNC processes receive contributions from tree level Z° and Higgs 

exchanges and bosons couple to right-handed quarks. In a particular version 

£T) ■ of this model in which the Yukawa matrix X D in the heavy down fermion sector 

is unitary, X u = 1 and M = Mvf is fixed, only three real and positive definite 

.Iph ' parameters describe New Physics (NP) contributions to all AF = 2 and AF = 1 

observables in K and B S ( i systems once the known quark masses and the CKM 
matrix are correctly reproduced. For M > 1 TeV NP contributions to B® d — B® d 
mixings are found to be very small. While in principle NP contributions to £k and 
AF = 1 processes could be large, the correlation between ek and Kl — > /U + 
eliminates basically NP contributions to ek and right-handed current contributions 
to AF = 1 FCNC observables. We find CMFV structure in B s d decays with 
B(B St d — > /*t + /U - ) uniquely enhanced for M = 3 TeV by at least 35% and almost up 
to a factor of two over their SM values. Also B(K + — > -k^vv) and B{Kl -4 iftvv) 
are uniquely enhanced by similar amount but the correlation between them differs 
from the CMFV one. We emphasize various correlations between K and B Sf d 
decays that could test this scenario. The model favours 7 ~ 68°, | V^f, | ~ 0.0037, 

S^Ks ~ °- 72 > %> ~ °- 04 and 42 x 10_9 < B ( B s "4 < 5.0 x 10- 9 for 

M = 3-4TeV. 



1 Introduction 



In a "minimal theory of fermion masses" (MTFM) [T] the fermionic content of the SM 
is extended by heavy vectorlike fermions, with flavor-anarchical Yukawa couplings, that 
mix with SM chiral quarks. Small SM Yukawa couplings arise then from small mixing 
angles between the heavy and light sectors. Although the hierarchical structure of the 
mixing is put in by hand, this model can be regarded as an effective description of the 
fermionic sector of a large class of existing flavour models and thus might serve as a 
useful reference frame for a further understanding of flavour hierarchies in the SM. 

This mechanism of quark mass generation implies corrections to the SM quark couplings 
to W , Z° and Higgs so that FCNC processes receive contributions from tree level 
Z° and Higgs exchanges and couple to right-handed quarks implying the presence 
of right-handed currents in a number of observables. The first question then arises 
whether this framework, while generating measured quark masses and CKM parameters 
is consistent with the constraints from electroweak precision observables, tree-level decays 
and in particular FCNC transitions for heavy fermion masses in the reach of the LHC. 
The second question is whether there is a particular pattern of deviations from SM 
predictions for FCNC processes in this model that could further test it in the flavour 
precision era. 

In [1] we have derived general formulae for the modified couplings of the SM heavy gauge 
bosons {W^, Z°, H°) to quarks and we have discussed the bounds on heavy fermion 
masses from present data finding globally consistency with the data for heavy fermion 
masses in a few TeV range. 

One of the central formulae in P was the leading order expression for the SM quark 
masses 

mZ = ve?eX)&, (X = U,D), (1.1) 



where 



Q,U,D 

v'lTn ( L2 ) 

with rrii describing the mixing between the light and heavy sectors and Mj standing for 
the heavy fermion masses (see Section 2). The matrices X U,D are the heavy Yukawa 
couplings which in [1] have been assumed to be anarchical 0(1) real numbers. While 
such an approach allowed a first look at the phenomenological implications of this model 
it did not allow a study of CP violation and a meaningful identification of correlations 
between various flavour observables. 

The goal of the present paper is to present a more general study with X U,D being first 
arbitrary complex matrices for which we will perform approximate but analytic diago- 
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nalization. Subsequently we will consider a specific structure for A ' requiring them 
to be unitary matrices. This will allow us to simplify the expressions of interest and 
significantly reduce the number of free parameters. In this manner a rather transparent 
phenomenology will follow. In particular we will find a number of correlations between 
different observables that can be tested in the future. Moreover it will turn out that 
in the simplest version of the unitary model (UM) in which X u = 1, to be termed the 
trivially UM (TUM), very definite predictions for particle-antiparticle mixing and the 
deviations from the SM expectations for rare K and B s ^ decays are found. 

Our paper is organized as follows. In Section |2] we summarize briefly those ingredients 
of the MTFM that we will need in our analysis. In particular we discuss the case of 
unitary heavy Yukawa matrices which specify the UM. In Section [3] we will perform 
the diagonalization of the induced quark mass matrices that will allow us very general 
expressions for quark masses, mixing angles and the couplings of SM heavy bosons to 
quarks in terms of the fundamental parameters of the model. We will then apply these 
formulae in the case of the UM and TUM in which we set X u = 1. The TUM has only 
three free real and positive definite parameters after all the CKM parameters and quark 
masses have been determined and a common mass M for heavy fermions has been fixed. 
This makes TUM very predictive. 

In Section H] we present the formulae for the effective Hamiltonians governing particle- 
antiparticle mixings K° — K° and B^ s — Bd,s that in addition to the usual SM box 
diagrams receive tree-level contributions from Z° and H° exchanges. Only induced tree- 
level Z° exchanges are relevant in MTFM. We also give a compendium of formulae 
relevant for numerical analysis of AF = 2 observables. In Section [5] the effective Hamil- 
tonians for s — > di/u, b — > qvv and b — > q£ + £~ (q = d, s) transitions are given. There 
we also comment on the contributions of W ± , Z° and H° to B — > decay which in 
the TUM for M of order few TeV turn out to be negligible. In Section [6] we calculate 
the most interesting rare decay branching ratios in the K and B meson systems, in- 
cluding those for the processes K + — > 7i + uu, Ki — > 7i vv, — > n°£ + £~ , B — > K^vv, 
B — > Xs dvv, B s d — > B + — > t + v t and Kl — > fi + {i~. Finally constraints from 

b — > s£ + £~ transitions like B — > K*(K)£ + £~ decays are discussed. 

In Section [7] we will search for a global structure of NP contributions in the TUM in 
order to understand better the numerical results presented in the subsequent section. 
This search is accompanied by the following questions: 

• Can the so-called Sk^S^k s tension present in the SM be removed in the TUM? 

• Which values of \V u b\ and of the angle 7 in the unitarity triangle are favoured 
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by the TUM on the basis of FCNC processes and what are the implications for 
B+ -> r+z/ r ? 

• What are the implications for the mixing induced CP-asymmetry and the 
branching ratios for B s ^ — > fi + [x~ and how this model faces new results from 
LHCb [2] ? 

• What are the implications for K + — > 7i + uu, Ki — > ifivv and Ki — > fi + fi~ decays 
that can probe high energy scales beyond the reach of the LHC [3]. 

In Section [8] a detailed numerical analysis is presented. In particular we study the 
correlations not only between various AF = 1 observables but also between AF = 1 
and AF = 2 observables. We summarize our results in Section [9j Few technicalities are 
relegated to appendices. 



2.1 Basic Ingredients 

The aim of this section is to briefly review the most important ingredients of MTFM. A 
detailed theoretical discussion is presented in pp. 

In this minimal model we add heavy vectorlike fermions that mix with chiral fermions. 
The hierarchical structure of SM fermion masses can then be explained through mass 
hierarchies entering the mixing pattern. We reduce the number of parameters such that 
it is still possible to reproduce the SM Yukawa couplings and that at the same time 
flavour violation is suppressed. In this way we can identify the minimal FCNC effects. 
The Higgs couples only to vectorlike but not to chiral fermions, so that SM Yukawas 
arise solely through mixing. However this mixing induces flavour violation already at the 
tree level proportional to v 2 /M 2 , where M is the mass of the vectorlike fermion. This is 
due to the fact that 

• SU(2)l doublets mix with SU(2)l singlets and 

• the Higgs couplings are no longer aligned with SM fermion masses. 

Now we specify the field content of the model focusing on the quark sector: we have 
three generations of chiral quarks 



2 The Model 




% = 



1,2,3 



(2.1) 
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and add for each of them a vectorlike pair of heavy quarks^ 



Um, U Li, D Ri , Du, Qb,i 



u Ri 
u Ri 



Qi 



1,2,3. 



(2.2) 



(2.3) 



The Lagrangian is then of the form (omitting kinetic terms) 

-C = h\%Q Li U Rj + h\°Q Li D Rj + M?U Ll U Rl + M?D Ll D Rl + M?Q Rl Q Ll 

+ mfU Li u Ri + mfD Li d Ri + mfQ Ri q Li + h.c. 

where we assumed that m x and M x (X = U, D, Q) can be diagonalized simultaneously. 
Then M x and m x are diagonal matrices with positive entries. Instead of m x we will 
use e x defined in fll.2p which are also real and positive. In order to get the low-energy 
effective Lagrangian that contains only SM fields we proceed in two steps: First we 
integrate out the heavy states using their equations of motions in unbroken SU{2) L x 
U(1)y- If we relax the assumption of aligned m x and M x we additionally have to 
diagonalize the mass matrix. As eigenstates we then get three heavy vectorlike quarks 
and three massless states which can be identified as SM quarks. Second we include 
EWSB and redefine light fields to get canonical kinetic terms. The relevant part of the 
effective Lagrangian in the mass eigenstate basis for gauge bosons and H° but in flavour 



basis for quarks reads [T 
9 



174 GeV) 



fi -/charged 



k- C + ^neutral 

H 



+ u Li m^u Rj + d Li m®d R j + -j= (u Li y^u R j + d Li y?d Rj ) + h.c. 
where the quark masses and Yukawa couplings are given as (summation over k 



-Q-X\X 



~ X \ X 

\k i 



III 



X 



{A x ) v e$e x \ x 



'ik k c j ^kj 

The charged and neutral current read 



( A X\ p Q F X \ X 

\ R/kj » k ik 



J charged 


= ULi 




J neutral 


= u Li 


$ij 




-d Li 


5 ij 



Yd^+umiA^^Ydnj 



.2 jfJ. 
wJclmag 



(2.4) 

1,2,3) 
(2.5a) 

(2.5b) 
(2.6a) 
(2.6b) 



*As for scales above the heavy quark masses the model contains 18 dynamical quarks, our model is 
not asymptotically free but the one-loop beta function at this scales is almost zero. We anticipate that 
the model being embedded in a larger gauge group will eventually be asymptotically free. 

2 The masses m^- D in £ c ff should not be confused with m U ' D 'Q in C, but as we traded the latter 
masses for e^ ,D,< ^ no problems of this sort should arise in what follows. 
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with (X = U,D) 



(^f)ij - MjfMg £i £j ' ~~ j^Q j^Q 8 ? £ f X i^k* , ( 2 - 7a ) 

(^n, = ^^X.A-, (2.7b) 

^ = -7T=Hnt> M fc x = M fc x (l + efef). (2.7c) 

We note that ^4^ ,£> enter both the charged and neutral currents. Explicit Feynman rules 
in mass eigenstate basis for both gauge bosons and fermions will be derived in Sec. [3j 



2.2 The Unitary Model (UM) 
2.2.1 Basic Assumptions 

We will now discuss an explicit model based on the Lagrangians (12. 3 p and (12.41) and the 
following two assumptions: 

• Universality of heavy masses: 

M® = Mf = Mf = M (2.8) 

for all i. The crucial assumption is that the matrices above are proportional to a 
unit matrix. The assumption that the masses in these three matrices are equal to 
each other is less important but helps in reducing the number of free parameters. 

• Unitarity of the Yukawa matrices X u and X D . While other forms, like symmetric 
or hermitian matrices could be considered, the unitarity of Yukawa matrix once 
assumed is valid in any basis and simplifies in a profound manner the charged and in 
particular neutral current interactions written in terms of fundamental parameters 
of the model. As we will see below, the resulting structure of interactions resembles 
the one in Randall- Sundrum (RS) scenarios. 



2.2.2 Parameter Counting and Explicit Parametrizations of A and A 



We have the following parameters: 

• 10 real parameters: M and e9' u ' D for % — 1, 2, 3. 
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A and A being unitary matrices have each 3 real parameters and 6 phases. 

Phase redefinitions of Ur, Dr and Ql allow, as clearly seen from the first two 
terms in the basic Lagrangian, to remove 8 phases. We are thus left with 



16 real parameters and 4 phases. 



(2.9) 



Subtracting six quark masses and four parameters of the CKM matrix our model has 
the following number of new parameters 



7 new real parameters and 3 new phases. 
The parameters in (12. 9p can be distributed as follows 



(2.10) 



10 real parameters: M and e 



Q,U,D 



for i = 1, 2, 3. 



As phase redefinitions of Ql can be made only once, X u and X D will have different 
number of phases but the same number (3) of real parameters. 



For X u and X D one can use explicit parametrizations of unitarity matrices known from 
the Standard Model (CKM matrix) and the corresponding matrix in the Littlest Higgs 
model with T-parity (LHT) describing the interactions of quarks with mirror quarks pE]. 

In the minimal version of MTFM we use for X D the CKM-like parametrization 



/ 



A 



D 



„d „d 
c 12 c 13 



„d „d 



d „d „d AS d „d „d 



~ S 12 C 23 C 12 S 23 S 13 e 



a d „d 
5 12 c 13 

n d „d d A5 d 



C 12 C 23 S 12 S 23 S 13 e 



„d e -iS A 

d „d 
6 23 c 13 



\ a d „d „d „d Q d p i8 d Q d „d Q d „d „d „i5 d ^d ^d 

\ ^12*23 L 12 c 23 6 13 e *23 c 12 & 12 L 23*13 e 



C 23 C 13 



(2.11) 



/ 



with cf- = cos 9f, and sf, = sin#£-. c£ and sf, can all be chosen to be positive and 5 d 
may vary in the range < 5 d < 27T. 

Following [1] the mixing matrix A^ can be conveniently parameterized, generalizing the 
usual CKM parameterization, as a product of three rotations, and introducing a complex 
phase in each of them, thus obtaining (c^- = cos 6^, = sin. 9^) 





(l 





x u = 





c 23 






— S2 3 e 4<5 23 



\ 



23 



/ 



C 



13 





1 

V-s? 3 e"5» 



'13 







( 


c 12 













Cl2 


Cl3 J 




K 





l) 








(2.12) 
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Performing the product one obtains the expression 

/ <-12 c 13 6 12 c 13 e 6 13 e > 

_S 12 C 23 e ~~ C 12 S 23 S 13 e C 12 C 23 ~~ S 12 S 23 S 13 e S 23 C 13 e 

\ S 12 S 23 e ~~ C 12 C 23 S 13 e _C 12 S 23 e ~~ S 12 C 23 S 13 e C 23 C 13 / 

(2.13) 

which completes the explicit parametrization of the model. 

However, even without looking at details the property of the unitarity of \ u and X D 
combined with the universality of the heavy fermion masses implies a very transparent 
structure of corrections to weak charged and neutral currents before one rotates to the 
mass eigenstates for SM quarks. 



2.2.3 Implications 

The general expressions for the gauge couplings simplify considerably. Indeed we find 



M 2 1 3 v w « M 2 



In order to obtain these expressions, without loosing the generality, we have modified 
the universality assumption in ( 12. 8 j) so that it is valid for and not Mj. As we will see 
in our phenomenological analysis Mj.,2 = Mi^ to an excellent accuracy but M3 7^ M3 
and this redefinition has to be made only in the case of M 3 . 

Inserting these expressions into neutral and charged currents we note that 

• The neutral current is diagonal in the flavour basis but the universality of inter- 
actions is broken by different values of £®' U,D for i = 1,2,3 that are necessary 
in order to explain the mass spectrum of quarks. Such universality breakdown is 
known from the RS scenarios where the universality is broken by the bulk masses. 
After rotation to mass eigenstates we will obtain FCNC currents at the tree level 
but they will be partly protected by the difference in masses. This structure is 
analogous to the RS-GIM [5]. 

• The left-handed charged currents are also diagonal within generations but the 
right-handed ones not. This will generally imply that the RH mixing matrices in 
the mass eigenstate basis will differ from LH-ones. 
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An interesting structure arises in this model also in the case of mass matrices and 
Yukawa couplings. Defining 



ij i 3 ij" 



(2.15) 



we find 



mE = v 



wPa = v 



A 



V 



ij ~ 2AP^ j ' ij 



2M 2 



2M 2 



^(f^\ 2 \ u — 



+ 0(v 5 /M A ) 
+ 0(v 5 /M*) 



llg 



vg 



(2.16a) 



(2.16b) 



Note that again if the ef' u,D where independent of the index % the Higgs couplings 
would be aligned with quark masses and we would not have tree level Higgs FCNC 
in the mass eigenbasis. 

The matrices \ u ' and A D ' in Eq. ( 12. 1 5[) are not unitary with the departure from 
unitarity given by the hierarchy of £j parameters. 



In order to simplify the notation in the rest of the paper we will suppress the bars and 
work with ef and M which really represent ef and M. In this context it is important 
to keep in mind that ef < 1 (see I2.7cl) so that with this notation in the subsequent 
formulae 

ef < 1. (2.17) 
2.2.4 Trivially Unitary Model (TUM) 

In our first phenomenological analysis it will be useful to consider a special case of the 
Unitary Model in which the Yukawa matrix X u is trivial: 

X u = l. (2.18) 

In this model, to be called TUM in what follows, we have then 

13 real parameters and 1 phase. (2-19) 
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This means that TUM has only four new parameters 

4 new real parameters and new phases. (2.20) 

Thus in the TUM CP violation is governed by a single CP phase S d in the Yukawa matrix 
\ D , which as we will show is equal to the CKM phase. Yet the phenomenology of this 
model will differ from the SM one due to new contributions from tree-level diagrams 
that contain new flavour violating parameters absent in the SM. In fact in this model it 
is very easy to express most of the parameters of the model in terms of the SM quark 
masses and the mixing angles of the CKM matrix. The remaining four free parameters 
can be chosen to be 

M, e%, 4, 4. (2.21) 

Note that all these parameters are positive definite and sf 3 and sf 3 are smaller than 
unity due to the unitarity of \ D . We will also see that in order to obtain the correct top 
quark mass for M of order few TeV one has 0.80 < s® < 1.0. Further restrictions on 
new parameters come from FCNC transitions. 

In order to perform phenomenology we have to rotate the quark fields to the mass 
eigenbasis and express ten fundamental parameters of the model in terms of six quark 
masses and four parameters of the CKM matrix. This we will do in the next Section. 

3 Diagonalization of the Quark Mass Matrices 
3.1 Preliminaries 

In this section we give general analytic expressions for the leading order rotation matrices 
V* R {X = U, D) that diagonalize the mass matrices = vX^efef via 

vFm x Vg = m$ as , X = U,D. (3.1) 

The CKM matrix is then given as 

Vckm = (Vf ) t V?. (3.2) 

With this at hand we can derive the flavour changing Z, W and H couplings in the mass 
eigenstate basis. 

In what follows we do not yet specify the parametrization of Yukawa matrices \ u and 
\ D so that the formulae in this section unless explicitly mentioned apply to MTFM at 
large. 
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Introducing 



e M 



e M 



M = Q,U,D 



(3.3) 



one can rewrite the mass matrices as 



ve q e x \ x e q e x 

V£ 3 £ 3 A ij £ i3 £ j3 ■ 



(3.4) 

The diagonalization is done as an expansion in powers of a small auxiliary parameter e 
that first can be set equal to the Cabibbo angle e = 0.23 §. To this end it is useful to 
have a rough estimate for the squared ratios ef- . We find 



-Q\2 



M \2 

-13J 



Kb ~ e 



mz 



Q\2 



T/2 ^4 



m 



c6 
2 1 



rsj rsj f 

ml VI 



( £ U)2 ~ e 



( -D \2 m d 1 2 
( £ 13) ~^2^2-~ e 



mz 



0.8-H1.6 



(3.5) 
(3.6) 

(3.7) 



We conclude that except for (e^) 2 , which appears too large to be considered as expansion 
parameter, we can expand in the remaining ratios. 

In our results for the mass eigenvalues, CKM matrix and rotation matrix we use the 
shorthand notation defined in Appendix |A] 

3.2 Quark Masses 



The mass eigenvalues in the up and down sector are given as 



m b 



m. 



ve$ef> 



1 x D 1 

QD \ A 22\ 



m d = ve^e 1 



det A 

I \D I 
l A 22l 



m r 



m,. 



Q u 

V £ 3 £ 3 



A33 = vefe%Kt 



ve?e° 



U \ A 22\ 



det X u \ 

\\U I 
l A 22l 



(3.8a) 
(3.8b) 

(3.8c) 



where for later convenience we have introduced parameters «j that collect the dependence 
on Yukawa couplings. This structure is very similar to the RS scenario (see 016] and 
references therein). The fermion shape functions f9' u,d correspond to our F 9 ,U ' D but we 
have here no exponential hierarchy but rather a powerlike one. Moreover, our calculation 
is a bit more general as we did not expand in e^. 



3 We emphasize that this value is indicative in order to perform perturbative expansion but has not 
been used in our numerical analysis. 
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3.3 Rotations to Mass Eigenbasis 



We parametrize the rotation matrices in the up sector in the following manner: 



/ 



i 



E?r,U 



12"1 



44\ 



-4(4)* 



£ 23 U 3 



\ £ 13 



4K) 



-£&(«?)• 



4«f 



-(7 /„,ii\* 



L 
3 

1 y 



£ 23 M 3 



/l \ 

e ifcc7 

\0 e icc/ 7 

/l \ 

e i6c7 

\0 e iCu J 



(3.9) 



(3.10) 



-eUniT 1 y 

with the 0(1) coefficients wf that are listed in Appendix [A] and phases that depend 
on the Zi defined below in (13 . 1 6[) according to 



o ib U 



Z l\ 



Z-i Zr> 



\ Z l\ \ Z 3\ 



(3.11) 



In the unitary model the coefficients u { ' can be written as functions of angles and 
phases. 

In the down sector rotations are parametrized as 



D 



D 
R 



/ 


1 


f q d L 




-4(4)* 


1 


V 


Q jL 

fc 13 U 4 


— 4(^3 




/ 1 






£ 12^3 


df 








ts ( 


jL,R 1 

i { (see 


Appendi 



c 13 u, 2 

^23^3 
1 



d R \ 



F D „ 



/l \ 

e if>D 
\0 e iCD 

/l \ 

e ibD 

ic D 



Jbr 



o lC U 



\0 & 

Dhases aci 

z\ z\ 



(3.12) 



(3.13) 



7 



z i\ 



\ Z l\ \ Z 3\ 



The translation to the notation used in [B] can be found in the Appendix 
our Vr is more general. 



(3.14) 
However 



3.4 CKM Matrix 

Using the rotation matrices derived above we can calculate the CKM matrix in (13. 2p . 
Expanding again in e and guided by ( 13.5l) -( |3~7l) we get 

S12C13 = A = \zt\ , S23C13 = AX 2 = \z 3 \ , si 3 e~ lS = AX 3 (p - ir]) = z 2 , (3.15) 



11 



with 



Z2 



fc 12 



'13 



"23 



X D 

A _12 
A 22 



A 12 
A 22 



4 (4 - «1 ) = 4"l2 



x D x^ 
\£> x^ 

A33 A 22 



X* 7 X 15 

A 12 A 23 
A 22 A 33, 



X D 
A 23 



^23 _ ^23 \ _ Q f,L _ L\ = 
x n \U I ~ £ 23 \ a 3 u 3j — 



A 33 



X 17 
A 33 



(3.16a) 



(3.16b) 



(3.16c) 



where for later convenience we have introduced parameters aij that collect the depen- 
dence on Yukawa couplings. 

Analogous formulae can also be found in [6] in Eq. (3.9), (3.10). The relation between our 
notation and the one used in [6] is summarized in Appendix From these three elements 
we can extract all four parameters of Vckm and then use the usual PDG convention 

/ C12C13 Si 2 ci3 SMe~ lS \ 



V ( 



CKM 



C12C13 S12C13 s 13 e 

S12C23 - Ci 2 S23Sl3e* 5 C12C23 - Si 2 S23Sl3e* 5 S23C13 
\ S 12S23 — c 12C23Sl3e i<5 — C 12 S23 — ■Sl2 c 23 s 13 ei<5 C23C13 J 



(3.17) 



The three parameters syi = \V us \/ci3, S13 = \V u b\ and S23 = |Vd,|/ci3 can be extracted 
from tree level decays and should not be sensitive to NP unless the charged currents are 
significantly modified by NP contributions. As we discuss in subsection 13.101 this turns 
out to be not the case in the TUM. 



3.5 Modified Z and W couplings 

While in the quark flavour basis the Z couplings were determined by the matrices A^ 1 ^, 
in the mass eigenstate basis they are governed by the matrices A* {X = U,D) given 

by m 

L = m cLg-^L m dLg > Aft = m ^ag-^R m di a g > (3.18) 

B? = V?B$Vg , Bl = V^B R V L X . (3.19) 

Here 

Bl = 4 diag (-£ , ^, ^) , Br — jjt diag ^) . (3.20) 

As seen in (I2.6ap and (I2.6bl) in quark flavour basis the same matrices A^ D appear in 
left-handed charged and neutral currents. However the rotation to mass eigenstate basis 
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affects charged and neutral currents in a different manner. The left handed W couplings 
are then governed by0 

Vckm , V^A U L VE = A u l Vcku , V^A°V° = VckmA? . (3.21) 

Induced right handed couplings of W are described in the quark flavour basis by A 1 ^ . 
In the mass eigenstate basis these couplings are given by 

A U R D = m U diag B U R D m^ , B™ = V^B R V L D . (3.22) 
As already stated in [1] the couplings scale according to 

^Iri^?- ( AX «lr^ (3 - 23) 

and for the right handed W coupling we get 

Inserting the rotation matrices V L ^ presented above and expanding in e we can calculate 
explicitly these flavour violating Z and W couplings. The analytic expressions for these 
couplings are rather complicated and can be found in Appendix iBl 

Now comes an important point. As seen in Appendix [B] in the MTFM all these flavour 
changing tree-level couplings (neutral and charged) depend on the same parameters that 
determine the SM quark masses and CKM parameters. However, generally the number 
of fundamental parameters is larger than ten and even after the correct quark masses and 
CKM parameters have been reproduced in this model, a number of free parameters will 
remain implying in principle non-MFV interactions. This point has been in particular 
stressed in a general context in [7]. We will investigate the size of these effects in TUM 
in Section El 

In Sec. 12.21 we showed that in the unitary model the gauge couplings - except A^ D for 
right handed W couplings - are diagonal in the flavour basis but different for the three 
generations. However, rotating to the mass eigenstate basis off-diagonal elements are 
generated: 

~Al = diag (e? , ef , ef) V? , (3.25) 



x 

2 



= j^V^ diag (ef \ e*\ sf ) V* . (3.26) 



4 The CKM matrix comes from the usual SM W couplings, but the second and third flavour violating 
coupling is new. 
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Explicit formulas are listed in Appendix [B] In order to simplify the notation we use in 
what follows shorthand notation: 

ef ee {ef)\ ef' ee (sf ) 2 , e? 4 ee ( e f )* (3.27) 



3.6 Final Results for A^ R (Z), A% R (W) 

In the expressions for effective Hamiltonians presented in the subsequent sections, the 
fundamental role is played by the couplings A l [ R (Z) that are defined as follows 

£fcnc(Z) ee [C L (Z) + C R (Z)} , (3.28) 

where 

C L (Z) = [Af{Z){s ltl P L d) + A b L d (Z)(bj,P L d) + A b °(Z)(b lfl P L s)} Z» , (3.29) 
C R (Z) = [Af{Z){- S1 ,P R d) + A h i{Z){h^P R d) + A b £(Z)(h»P R s)] Z" , (3.30) 
and A l [ R {Z) are the elements of the matrices Al jR (Z). 

The modifications of the W couplings to the SM quarks, A l l R (W), can be similarly 
summarized by the Lagrangian: 

A£(W) ee [C L (W) + C R (W)} , (3.31) 

where 

C L (W) = [A td (W)(i ltl P L d) + A£(W0(t 7 „P L s) + A^iW)^^)] W , (3.32) 

C R (W) = [A%{W)(t 7lt P R <t) + A^W^Prs) + A tb (W)(i^P R b)] , (3.33) 
and A l [ R (W) are the elements of the matrices Al^ r {W). 

In the mass eigenstate basis these couplings are given in terms of A* R and A V R D as 
follows 

A^(Z) = --^(i£%, (3.34) 

Al(Z) = ^-{A^ . (3.35) 



Al{W) 



2V2 



AYVckm ) + ( VckmA d 



ij V / 



(3.36) 



Al(W) = (A U R D ) ^. (3.37) 



^2 ; 



Note that 



Aiy^) = (^UZ)T, Al R {W) = (A i l R (W)T . (3.38) 
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3.7 Results for A l [ R {H) 

For completeness we list the flavour violating H couplings even if in the minimal version 
of our model the Higgs contributions to all processes considered by us are as expected 
negligible [I]. 

In analogy with the flavour violating couplings of Z Q we introduce 

£fcnc(#) = [C L (H) + C R (H)} , (3.39) 

where 

C L (H) = [A s L d (H)(sP L d) + A b L d (H)(bP L d) + A b L s (H)(bP L s)] H, (3.40) 
C R (H) = [A sd (H){sP R d) + A"(H){bP R d) + A b £(H)(bP R s)] H. (3.41) 

The Higgs couplings to light down-type quarks are given as = d, s, b): 

mf(l£)*,+ (i£)>?) , (3-42) 
(Afiymf + mf(ig)«) . (3.43) 

(3.44) 

We observe that relative to flavour violating Z° couplings, these couplings are suppressed 
for all SM quarks but top quark by roughly rrii/v where m 8 denote the masses of SM 
quarks. Consequently H° contributions do not play any role in the processes considered 
by us. 

3.8 Expressing Model Parameters in terms of Quark Masses 
and CKM Parameters 

For our numerical analysis it is useful to express analytically as much as possible the 
fundamental model parameters in terms of the measured quark masses and CKM param- 
eters so that the number of free parameters relevant for the analysis of various decays 
will be significantly reduced. Simultaneously correct values of quark masses and CKM 
parameters will be incorporated in our analysis automatically. We consider first the 
general case and subsequently show the case of the TUM where the analysis becomes 
very simple and transparent. 



Af{H) 
A%(H) 



V2v 
1 

V2v 
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3.8.1 General Case 



First from (I3.8P we can determine ef' in terms of quark masses, ef and the parameters 
in the Yukawa matrices that are hidden in Kj. We find 

D _ I m h v _ 1 m t 

sfvn t 

1 m r . 



1 


m b 


e Q 

fc 3 




1 


m s 


e Q 
-2 


VK S 


1 


m d 



-=^— • 4 = 5^' (3 - 45b) 

e? = ^ , e u 1= ^. (3.45c) 

ef vK d £ f vk u 

Next using (13.161) we can determine ef and ef in terms of S13, S23, ef and the parameters 
in the Yukawa matrices that are hidden in aij. We find 

e? = e?T^, ef=ef^ (3.46) 

|«13| " F23| 

Finally writing 

ay = |ay|e"^ (3.47) 

we find the following two conditions 

2 , 013 = <$CKM = 1, (3.48) 

■S23 S 12C 13 





a 


13 




«23 


a 


12 



with 7 being one angle of the unitarity triangle that up to the sign equals the phase of 
V ub . 

In this manner we could directly express analytically the fundamental parameters of the 
model in terms of six quark masses and three mixing angles that can be determined 
from high energy collider experiments and tree-level decays, respectively. With two 
conditions in (I3.48P we can fix one phase and one mixing parameter in the Yukawa 
matrices in terms of the remaining five mixing parameters in these matrices and three 
phases. The additional two real free parameters are ef and the common mass M of the 
heavy fermions. 

At this stage one comment should be made. This procedure can be straightforwardly 
executed numerically if one is allowed to set = in the expression for A^, which 
appears to be a good approximation if the diagonal terms in A D are largest. Otherwise 
these formulae could be used iteratively setting first = 0. 
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Figure 1: Correlation between s| 3 and sf 3 in the TUM for \V u b\ = 0.0037 for different 
sf 2 : 0.05 (dashed red), 0.2 (dotted blue), 0.4 (dot dashed green), 0.6 (solid purple), 0.8 
(solid cyan), 0.95 (solid magenta). 



3.8.2 Trivially Unitary Model 

In this case the number of free parameters is decreased to four that are listed in (I2.2ip . 
Even more importantly the two conditions in (I3.48j) simplify considerably 



In this manner all conditions from quark masses and the measured CKM matrix from 
tree-level decays can be automatically satisfied and once these parameters have been fixed 
all the observables analyzed by us are given entirely in terms of the four real parameters 
listed in (I2.2ip . In Fig. [T]we show the correlation between sf 3 , s 23 and sf 2 implied by 
the condition in (I3.49[) using \V u b\ = 0.0037. As we will include in our analysis many 
observables we expect a number of correlations between them in this scenario. We also 
emphasize that in this model CP-violation is governed entirely by the CKM phase but 
as in addition to CKM parameters new flavour parameters are present we will find some 
deviations from the correlations present in the SM. 



c d \„iS n d -t-d I a d Q d 
6 23 | e c 23 t 12 ' 6 23 6 13 



5 = 5 C KM = 7- 



■S23 S 12C 13 




(3.49) 
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Moreover in TUM 1 and we find in particular an important relation 

s^ = ^Ml. (3.50) 

For our nominal value M = 3TeV one finds m t (M) = 139 GeV. Taking into account 
that £-§ < 1 and v = 174 GeV we find the allowed range for 

0.80 < ef < 1.0. (3.51) 

As this parameter plays a prominent role in our phenomenology, this range, following 
from the desire to obtain correct top-quark mass, has significant impact on possible devi- 
ations of FCNC observables from SM expectations. As m t (M) decreases with increasing 
M, the lower limit in f !3.5ip also decreases with increasing M but this decrease is only 
logarithmic and in the range of M in the reach of the LHC, this dependence can be 
neglected. 

The beauty of this simple scenario lies in the fact that at the end, when the data improve 
and the theoretical uncertainties in FCNC observables will be reduced, we will be able 
to fix completely the Yukawa interactions in the heavy vectorial quark sector and the 
mixing parameters e9< u ' d that connect this sector with SM quarks. 

3.9 Couplings A* R in the TUM 

The general expressions for A* R and A V R D even in the unitary model are rather lengthy 
and have been presented in Appendix [B] However, they are simpler in the TUM and we 
list them here. 

In the TUM the up- type couplings simplify considerably, since = and K u ,c,t = 1- 
Consequently we get diagonal up-type couplings: 

% = ^iag ( % % ™k) = ^Uag (ef , ef , ef) , (3.52) 



A « = ii^ diag ! S' ^) = w diag £ " £ ^ ■ (3 - 53) 

1V1 Q V^l £ 2 £ 3 / 1V1 Q 

This means that in the TUM there are no-tree level contributions to FCNC processes in 
the up-quark sector, e.g. in D° — D° system. If one day NP effects in this system will 
be convincingly shown to be present, the condition X u = 1 will have to be removed. 

The down-type couplings are non-diagonal and given as 



A ° L = ^ckmcM? (e? 2 , ef, ef) V C km , (3.54) 

= TH m Sag^CKMdi a g ( ~Q2> ~Q2' ~Q2 ) ^CKM^diag > (3.55) 
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These formulae are very transparent. Indeed all phase matrices drop out. This means 
that all information about FCNCs in the TUM are encoded in the ef. We again observe 
that for 

e? = e? = e? (3.56) 

the coupling matrices A® and are diagonal and there are no FCNC transitions 
mediated by Z° at the tree level. However, the hierarchical structure of the CKM matrix 
breaks this relation as in MTFM 

£ Q £ Q £ Q 

—q = |Kis|^12) —q — \Vub\X13, —q = \V c b\X 2 3, (3.57) 

e 2 £ 3 £ 3 

where we have introduced quantities that will be useful later on. In [1] Xy have been 
assumed to be close to unity. In the case of arbitrary X D and X u they can differ signifi- 
cantly from unity and from each other. But when all constraints are taken into account 
the equality (13.561) remains badly broken and tree-level FCNC transitions mediated by 
Z° are present. 



3.10 Modification of W ± Couplings in the TUM 

We analyse next the modification of the charged couplings in the TUM. Using the general 
formula for the corrections to left-handed couplings in ( 13.361) we find the effective CKM 
matrix 

[iTL = C'(i4 £f,) ' (3 - 58) 

The effective CKM matrix is clearly non-unitary, a property know from other studies 
in which heavy fermions mix with the SM quarks [5HTT]. See for instance [H], where a 
detailed study of this effect in the context of a RS scenario has been presented. With 
v = 174 GeV, M = 3.0 TeV and ef < 1 these corrections are at most at the level of 
0.5% and thus negligible. Therefore, in our numerical analysis it is legitimate to use 
the unitarity of the CKM matrix and neglect the corrections to left-handed couplings of 
W ± . 

For the right-handed W ± couplings we get 

^R D = TTJ^diagdiag ( -gj, -55, -gj J V C KMm^ g (3.59) 
and consequently using (13.371) 



A « (,v) = -jfW^?-' (3 - 60) 



19 



where mf and are SM quark masses normalized at M. 

Using this formula we can check, whether the presence of right-handed currents could 
help in explaining the differences between the determination of \V u b\ in exclusive semi- 
leptonic decays, inclusive B decays and B + — > t + u t as proposed in [12j[T3] and analysed 
in detail in [MfT5] . 

Following [H] we find for exclusive semileptonic decays the effective V u b 

KT cl = C + V ub , (3.61) 

while in the case of B + — > t + u t 

V: b = C.V ub , (3.62) 

where 

C7 ± = l±^»-p. (3.63) 

We find that the sign of the correction in C + is opposite to the one required for the 
explanation of the difference between inclusive and exclusive \V u b\ determination but the 
correction turns out to be so small that it can be neglected. Similar comment applies to 
C_. Indeed using (13 . 5 7f) we have for M = 3.0 TeV 

C ± = 1 ± m :™ b m 1 « 1 ± 2^^10- 5 (3.64) 

M2 (^i 2 3 lK fe | 2 ) efXl, 

However, our analysis of FCNC processes will imply X 13 > 2 rendering these effects to be 

totally negligible. Recent decrease of the experimental branching ratio for B + — > t + v t 

[TBlfTT] makes the case for RH couplings of much weaker than two years ago anyway. 

With all this information at hand we are now ready to turn our attention to FCNC 
transitions in this model. We should emphasize that the formulae given in the next two 
sections apply to MTFM at large and only when the couplings A l [ R (Z) are specified to 
the TUM, one obtains the formulae specific to this model. 



4 AF = 2 Transitions 



4.1 Standard Model Results 



The dominant contributions to the off-diagonal elements M\ 2 in the neutral Kaxid B q 
meson mass matrices come from SM box-diagrams. They are given as follows 



(MS 



SM 



12^2 



F^B K m K M^ [A^SqOc) + A^ 2 S , (^) + 2\ c \ t i l3 S (x c , x t )] ,(4.1) 



127SM 



12tt 



-F 2 B B Bd m Bd M, 



AS 9) ) VBSo(x t ) 



(4.2) 
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where Xi = m^/M^ and 

A, = V*V id , \f = V*V iq (4.3) 

with Vij being the elements of the CKM matrix. Here, So(xi) and Sq(x c , Xt) are one-loop 
box functions for which explicit expressions are given e. g. in J4]. The factors rji are QCD 
corrections evaluated at the NLO level in [18H22J. For 771 and 773 also NNLO corrections 
are known [231124"] . Finally Bk and B Rq are the well-known non-perturbative factors. 

It should be emphasized that in the SM only a single operator 

QX hh (K) = {s llx P L d) (s^P L d) (4.4) 

and 

Q^ LL (B q ) = {h,P L q) (h^Lq) (4.5) 
contributes to and M\ 2 {q = d,s), respectively. Moreover complex phases are only 
present in the CKM factors. 

Our next goal is to generalize these formulae to include the new tree level contributions 
from Z exchanges as shown in Fig. [2j We will see that three distinct new features will 
characterize these new contributions: 

1. The flavour structure will differ from the CKM one. 

2. FCNC transitions will appear already at the tree level as opposed to the one-loop 
SM contributions in (gUD and ffl~2]) . 

3. In addition to QY LL (K) and QY LL (B q ) (with q = d,s) new operators will be present 
in the effective Hamiltonians in question. 



4.2 Tree Level Z Contributions 

We begin our discussion with the tree level Z exchanges contributing to AS = 2 transi- 
tions in Fig. [2j Analogous diagrams contribute to B® s — B® s mixings. 

The FCNC Lagrangian is given in fl3.28p - fl3.30p in terms of A l [ R (Z). These couplings 
are complex quantities and introduce new flavour and CP-violating interactions that can 
have a pattern very different from the CKM one. Explicit expressions for A l [ R (Z) in 
MTFM are given in Section 13.61 and Appendix [B] 

The diagrams in Fig. [2] lead to the following effective Hamiltonian for AS* = 2 transitions 
mediated by Z: 

Mff =2 ]z = ^72 \H d (Z)} 2 (s^P L d) {JrfPiA + [A#(Z)] 2 (s lfl P R d) (sYPRd) 



2M Z . 



+ 2Af(Z)A s «(Z) (s^P L d) (sYPnd) . (4.6) 
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Figure 2: Tree-level flavour changing Z° contribution to K — K mixing (the diagram 
rotated with 90° also exists). 

For the B® s — B® s Hamiltonians one has to replace u sd" by "fed" and "fes" , respectively. 

The Hamiltonian in ( 14. 6 j) is valid at scales O(Mz) and has to be evolved to low energy 
scales \x = (9(3 GeV), \x = 0(mb) at which the hadronic matrix elements of the operators 
in question can be evaluated by lattice methods. The relevant anomalous dimension 
matrices necessary for this renormalization group evolution have been calculated at two- 
loop level in [25l[26] and analytic formulae for the relevant QCD factors analogous to r\i 
in (14. ip and (14. 2 ft can be found in [2?]. The latter formulae do not include the 0(a s ) 
corrections to the Wilson coefficients of the relevant new operators at fiz = 0(Mj). 
Recently such corrections to AF = 2 observables originating in tree level contributions 
of gauge bosons and scalars have been calculated in [28] and we will include them in our 
analysis. 



In order to use the renormalization group analysis of [271128] we recall the operator basis 
used there: 

Qi LL = (s 7li P L d) (sYP L d) , (4.7a) 
Q7 RR = (s^P R d) (s^P R d) , (4.7b) 
Q\ K = {s llx P L d) (sYPnd) , (4.7c) 
Q\ K = (sP L d) (sP R d) , (4.7d) 

where we suppressed colour indices as they are summed up in each factor. For instance 
s^^Pid stands for s a7ll PLd a and similarly for other factors. 

A straightforward calculation gives us the effective Hamiltonian for AS = 2 transitions 
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in the basis (14.71) with the Wilson coefficients corresponding to = O(Mz) [28] 



(Af^))* V LW.. ^VLL , (A^(^)) 2 



^ eff Jz - 2M2 Cl {fIz)Ql + 2Ml Cl { ^ z)Ql 

Z 



VRR 



(4.t 



where including NLO QCD corrections [28] 



a. / „, Mi 11 



^(aO = C7™(A*) = 1 + ^ ^-2 log ^ + f ) • (1.9) 



^(,)-l + ^^lo g -/-g), (4.10) 



2 



C 2 L V) = ^(-61og^-lj. (4.11) 

The latter coefficients are also valid for B® — B® and B® — B® systems, where in the rest 
of the formulae sd should be replaced by bd and 6s, respectively. 

We should remark that also tree-level Higgs contributions to AF = 2 processes are 
present. However, as we have seen the corresponding flavour violating Higgs couplings 
to light quarks are very strongly suppressed and these contributions are totally negligible. 
We will not discuss them further. The general structure of such contributions including 
NLO QCD corrections can be found in [28] . 



4.3 Hadronic Matrix Elements 



In order to complete the analysis of AF = 2 processes we have to include renormalization 
group QCD evolution from the high energy scales down to scales at which hadronic matrix 
elements are evaluated. In what follows we will summarize the expressions for tree-level 
contributions to the off-diagonal elements M 12 for the AS* = 2 transition. Analogous 
expressions for B s d systems can be easily obtained in the same manner. 

In presenting our results we will use the so-called Pf QCD factors of [27] that include both 
hadronic matrix elements of contributing operators and renormalization group evolution 
from high energy to low energy scales. These factors depend on the system considered, on 
the high energy matching scale and on the renormalization scheme used to renormalize 
the operators. The formulae for these factors have been given in [27] in the MS-NDR 
renormalization scheme. This scheme dependence is canceled by the non-logarithmic 
0{a s ) corrections calculated in [28J and given in the previous subsection. The logarithmic 
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corrections calculated also there cancel the scale dependence of P? as demonstrated in 
that paper. 

The formulae for various contributions to M\i defined in the case of AS = 2 through 

(M*Y = (K°\n* s = 2 \K°) , (4.12) 

are easily obtained from the Hamiltonians presented above by replacing the operators 
by their hadronic matrix elements 

<<2?(Mz)> = ^j^PtM (4-13) 

with analogous definition for the B s d systems. We list now the resulting NLO expres- 
sions. For the Z contribution we find 128 



d d z (4 - 14) 

Z 

The relevant Wilson coefficients for AS = 2 are given in Eq. (14.91) with obvious replace- 
ments for B s l i systems. 

In order to find numerical values of P" and consequently (QKnz)) one needs the values 
of the corresponding non-perturbative parameters Bf defined in [27J. These are given in 
terms of the parameters Bi used in [29H3T] as follows: 



BY Lh M = B™-Qm>) = BM , (4.15a) 
£^ R (/io) = B B (no) , (4.15b) 
£? 2 LI Vo) = B 4 (fi ) . (4.15c) 

In the case of the K° — K° system, the values for Bi in the MS-NDR scheme have been 
provided in [29],[30] . We have just used the average of the results in [29l[30] that are 
consistent with each other. On the other hand the most recent results for the B® s — B® s 
systems [31] in the same renormalization scheme are given not for Bi but for B^ S F^ as 
this reduces the errors in (Q^(fi H )). In Tabled] we collect the values of 

Df(K) = Bf(fi L )F* for i = 4, 5 , (4.16) 
Df(B q ) = Bf(ji b )F* q for i = 4, 5 , q = d,s, (4.17) 

at the relevant scale /j,q given in the last column. 
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The effective parameters Bf(fii) are defined in the case of K° — K° mixing (//£ = 3 GeV) 
by 

3 eff/ \ — ( m K 



Bi{ni 



(4.18) 



m s y + m d (ii L ) 

with the most recent values for i? 4 5 at \i L = 3 GeV in the MS-NDR scheme given 
by [291130] 

B A = 0.76 ± 0.07, B 5 = 0.56 ± 0.06 . (4.19) 

In the case of B^ s — B® s mixings one has to make the replacements /i^ — > fit, = 4.2 GeV 
and change appropriately flavours 

2 

Bib*) (4.20) 



m B 



m b (fj, b ) + m d {n h ) 

with an analogous formula for the B° s — B Q S system. The values of weak decay constants 
and of meson masses required to obtain the hadronic matrix elements by means of (14.131) 
in the case of K° — K° system and analogous expressions for B® s — B^ s systems are 
given in Table HI For the SM contributions we use the necessary input from the latter 
table. 





ncff 
-1^4 


Df 




K°-K° 


0.675 


0.517 


3.0 GeV 


B° d -m 


0.093 


0.127 


4.2 GeV 


B° s -B° 


0.135 


0.178 


4.2 GeV 



Table 1: Central values of the parameters Di(M) in units of GeV 2 in the MS-NDR 
scheme based on f2^ l3U f for K° — K° system and 137]/ for B^ s — B^ s systems. The scale 
fio at which Di(M) are evaluated is given in the last column. 

Finally, we collect the values of {Qi(f^z)) contributing to (M 12 )z for Hz = ™>t i n Table [2J 
As the SM Wilson coefficients in the SM are also evaluated at m t we have also chosen 
this scale for the Wilson coefficients of new operators. As we include NLO corrections 
to the Wilson coefficients the same final results for the mixing amplitudes up to higher 
order corrections would be obtained if we used \iz = M z or any value 0(M w ,m t ). 

Concerning the new contributions to the Wilson coefficients of Q^ LL and the contri- 
butions from QY RR operators they can be included effectively by replacing the flavour 
independent S (x t ) in the SM formulae by the functions S(M) (M = K, B d , B s ): 

S{M) = S (xt) + [A5(M)]vll + [A5(M)] V rr. (4.21) 
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Wl \" b t!l 




K°-K° 


-0.11 


0.18 




-0.21 


0.27 


B° s -B° s 


-0.30 


0.40 



Table 2: Hadronic matrix elements (Q1{Tn t )) in units of GeV 3 at m t = 163 GeV. 



One finds |3j 
[AS(B q )] Vhh 



-\ 2 



4f 



[AS(AT)] V ll 



Af{Z) 



2 



A 



(K) 



4r 



(4.22) 



where gsM is defined in (15. 2 p and f = 1.068. [A5'(M)]vrr is then found from the formulae 
above by simply replacing L by R. The important new property is the flavour dependence 
in these functions and the fact that they carry new complex phases. 



4.4 Combining SM and Tree Contributions 

The final results for Af£, M? 2 and Mf 2 , that govern the analysis of AF = 2 transitions 
in the MTFM in question, are then given by 

M{ 2 = (M{ 2 ) SM + (M{ 2 ) z = (M{ 2 ) SM + (M{ 2 ) Np , {i = K,d,s), (4.23) 

with (Mi 2 ) SM given in (filjl -dO) and (M* 2 ) z in (Ol) . 



4.5 Basic formulae for A.F = 2 observables 

Having the mixing amplitudes M\ 2 at hand we can calculate all relevant AF = 2 ob- 
servables. To this end we collect below those formulae that we used in our numerical 
analysis. 

The Kl — Kg mass difference is given by 

AM K = 2 [Re {Mg) m + Re {Mg) m ] , (4.24) 
and the CP-violating parameter Ek by 



cxp 

26 



£K = V2(AM K ) CXP [lm ^ SM + Im ^ Np] ' (4 ' 25) 



where <p s = (43.51 ± 0.05)° and k £ = 0.94 ± 0.02 p|!53] takes into account that <p e ^ f 
and includes long distance effects in Im(Ti2) and Im(Mi2). 

For the mass differences in the B® s — B^ s systems we have 

AM, = 2 |(Af? 2 ) SM + (Mf 2 ) Np | (q = d,s). (4.26) 
Let us then write [31] 

M q 12 = (M« 2 ) SM + (M* 2 ) NP = (Mf 2 ) SM C^e***. , (4.27) 

where 

«)sM=l«)sMl e ^' ^ 22 °> ( 4 - 28 ) 

(M 1 S 2 ) SM = |(M 1 S 2 ) SM | e 2 ^ , p, ~ -1° . (4.29) 

Here the phases /3 and (3 S are defined through 

V td =\V td \e- 4f) and ^ = -|y ts | e -^ . (4.30) 

We find then 

AM, = (AM 9 ) SM C Bg , (4.31) 

and 

S^ Ka = sin(2/3 + 2(p Bd ) , (4.32) 

= sin(2|/3 s | -2^ B J, (4.33) 

with the latter two observables being the coefficients of sin(AMrft) and sin(AM s t) in the 
time dependent asymmetries in B^ — y tpKs and B° s — y ip<f), respectively. At this stage a 
few comments on the assumptions leading to expressions (14.321) and f!4.33|) are in order. 
These simple formulae follow only if there are no weak phases in the decay amplitudes 
for B® —y ipK$ and B° s —y ipcf) as is the case in the SM and also in the LHT model, where 
due to T-parity there are no new contributions to decay amplitudes at tree level so that 
these amplitudes are dominated by SM contributions [4j. In the model discussed in the 
present paper new contributions to decay amplitudes with non-vanishing weak phases 
are in principle present at tree level. However, as we demonstrated previously in the case 
of TUM these contribution can be totally neglected when calculating S^k s and S^. 

5 Effective Hamiltonians for AF = 1 Decays 
5.1 Preliminaries 

The goal of the present section is to give formulae for the effective Hamiltonians relevant 
for rare K and B decays that in addition to SM one-loop contributions could generally 
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include tree level contributions from the SM Z gauge boson, the tree level neutral Higgs 
H° and the corrections to SM one-loop contributions due to the modification of the W ± 
couplings. In the absence of tan j3 enhancement present in supersymmetric models and 
generally multi-Higgs models, the H° contributions can be neglected as their flavour 
violating couplings to quarks are strongly suppressed in our model and the couplings 
to leptons are strongly suppressed by small lepton masses. Similarly we find that in 
TUM the corrections from modified W ± couplings in loop diagrams governing in the SM 
rare K and B decays are much smaller than the tree-level Z° amplitudes. This can be 
expected on the basis of our estimates of such corrections in subsection 13.101 Therefore 
in what follows we will present only the effective Hamiltonians based on the SM loop 
contributions and the induced tree level Z° exchanges. 

The case of radiative decay B — > X s 7 is special as here there are no contributions in 
the SM and MTFM at the tree level and one has to check whether the modifications of 
W couplings and Z° exchanges in loop diagrams absent in the SM can have a visible 
impact on the branching ratio. In particular the presence of right-handed W couplings 
in this decay can lead to enhancement factors m t /mf,. Moreover as now W ± 1 Z° and H° 
couplings between light and heavy quarks are involved one has to check whether heavy 
quark contributions to B — >■ Xgj are significant. Following formalism for W ± and Z° 
contributions with both left-handed and right-handed couplings developed in [35] in the 
context of gauged flavour models and in [36] in the context of SU{2) L x SU{2) R x U{1) 
model, we have verified that the NP contributions from W and Z° exchanges to B — > 
X s 7 even in the presence of heavy quarks are negligible in the TUM. Extending this 
formalism to H° internal exchanges we have found that also these contributions are 
very small. As this analysis, including QCD corrections, is rather involved and the NP 
corrections in the TUM are negligible anyway, we will present it elsewhere in the context 
of a more general version of the MTFM. 

In [37] general formulae for effective Hamiltonians resulting from tree level neutral gauge 
boson exchanges with arbitrary masses and arbitrary left-handed and right-handed cou- 
plings have been presented in the context of the RS scenario. We could in principle apply 
them directly to our case. However, as in our case we have only one neutral gauge boson 
instead of three, we can write the relevant formulae in a simpler form, as done already 
in [3], than it was done in [37]. 
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Figure 3: Tree level contributions of Z to the s — > dvv effective Hamiltonian. 



5.2 Effective Hamiltonian for s — > dvv 

The effective Hamiltonian for s — > dvv transitions resulting from Z-penguin and box 
diagrams is given in the SM as follows 

Msm=4 E [KX^ t (x c ) + X t X(x t )](sj,P L d)(P a ,PL^)+h.c., (5.1) 

l=e,n,T 

where Xi = mj/Myy, K = Vi* s Vid an d V%j are the elements of the CKM matrix. X^ h (x c ) 
and X(x t ) comprise internal charm and top quark contributions, respectively. They are 
known to high accuracy including QCD corrections [381 - HO] and electroweak corrections 
|41[l4"2]. For convenience we have introduced 

V2 27r sin 6 W 

In the MTFM (15.11) is modified by tree-level diagrams in Fig. [31 The FCNC Zsd vertex 
has been given in (I3.28I) - (I3.30I) with explicit expressions for A| d (Z) and A S ^(Z) given in 
Section 13.61 and in Appendix iBl 

For the Zvv coupling we analogously write 

C VP {Z) = A?(Z)(9 lfl P L v)Z», A^(Z) = (5.3) 

A straightforward calculation of the diagram in Fig. [3] results in a new contribution to 

= -JfT 1 [&l d (Z)(s^P L d) + Af{Z){sYP R d)] (yi»P L v) + h.c. . (5.4) 

z 

Combining this contribution with the SM contribution in (15.11) . 
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we find 



[KI] K =9 2 sm E [KX^{x c )+\ t X L {K)]{s^P L d){9a^Pm) 

£=e,/i,T 

+ 9 2 sm J2 [^ K) X R {K)]{srP R d){9 af ,PLUt) + h.c. . 



(5.6) 



£=e,fi,T 

Here we have introduced the functions X L (K) and X R (K) 

XL{K) = x{Xt) + ^MlVZVu^ 

9 2 suM 2 z V t * s V td • 



X R {K) 



(5.7) 



(5.8) 



Finally the SM loop function X(x t ), resulting from Z-penguin and box diagrams is given 
as follows 



X(x t ) =T] X — 



x t + 2 3x t — 6 

+ t — In x t 



^ _x t — l (x t - l) 2 
r)x is QCD correction to these diagrams [381113] when m t = m t {m t ). 

Some comments are in order: 



T] X = 0.994 



(5.9) 



• In the SM only a single operator (sj^PLd) {y^^PiM) is present. This is due to the 
purely left-handed structure of SU (2) l gauge couplings. 

• In the MTFM in question also the operator {s^^ 1 PRd){i"-) ll PLi l ) is present, as Ar(Z) 
coupling matrices have non-diagonal entries. 

• On the other hand we will keep the neutrino couplings SM-like, that is purely 
left-handed. 

• As all NP contributions have been collected in the term proportional to , 
X^, NL (x c ) contains only the SM contributions that are known including QCD cor- 
rections at the NNLO level [39"P0"]. 

5.3 Effective Hamiltonian for b — > dvv and b — > svv 

Let us now generalize the result obtained in the previous section to the case of b — > dvv 
and h — > svv transitions. Basically only two steps have to be performed: 

1. All flavour indices have to be adjusted appropriately. 

2. The charm quark contribution can be safely neglected in B physics. 
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The effective Hamiltonian for b — > qvv (q = d, s) is then given as follows: 



IK 



VV]Bq 

effJ 



9 2 su E KV tb X L {B q )]{qrP L b){u a ,P L ^) 

£=e,fx,T 

+9 2 su E KV tb X R (B q )](q^P R b)(9 a »Pm) + h.c. 



with 



(5.10) 



X L (B q ) = X(x t ) 



Ag(Z) 



(Z) 



(5.11) 



Again all relevant A b £ R entries in the MTFM can be found in Section l3T6l and Appendix iBl 

Note that the functions X L>R (K) and X L ^ R (B q ) presented above depend on the quark 
flavours involved, through the flavour indices in the A l [ R (Z) = s, d, b) couplings and 
through the CKM elements that have been factored out. While in principle A % [ R (Z) 
could be aligned with the corresponding CKM factors, this is generally not the case and 
the functions in question become complex quantities that are flavour dependent. This 
should be contrasted with the case of the SM and CMFV models where the decays in 
question in the K, B^ and B s systems are governed by a flavour-universal loop function 
X(xt) and the only flavour dependence enters through the CKM factors. Consequently 
certain SM-relations and more generally CMFV-relations can be violated in MTFM. 
However, as we will see below in TUM they are satisfied with high precision in the B s d 
systems once all existing constraints on FCNC processes are taken into account. The 
case of K + — > ir + v& is different. 



5.4 Effective Hamiltonian for b — > dl + l and b — > s£ + £ 

The effective Hamiltonian for b — > q£ + £~ (q = d, s) can straightforwardly be obtained 
following the derivation of the effective Hamiltonian for s — > d£ + £~ transition, presented 
in [37] and properly adjusting all flavour indices. In addition, in contrast to the s — > 
d£ + £~ transition, now also the dipole operator contributions mediating the decay b — > 57 
become relevant. Explicit formulae for these contributions will be presented below. In 
the following we will denote the total contribution of the dipole operators to the effective 
Hamiltonian in question simply by "H e fr(& — > S7). As already mentioned previously in 
the TUM, this Hamiltonian is governed by SM contributions and its explicit form will 
not be given here. 
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The relevant Feynman diagram for the Z contribution, shown in Fig. HJ contains on the 
1. h. s. the vertex, which we already encountered in the case of the b — > svv decay The 
relevant FCNC Lagrangian for Zhs couplings has been given in f l3.28p - fl3.30p . For the 
£ + £~ vertex we write in analogy to (15. 3 p 

C d (Z) = [A^(Z)(£ lfl P L £) + A%(Z)(£ lfM P R £)] Z» , (5.12) 

where with definitions in (I5.2ip 

A^(Z) = -^-, A^(Z) = -^-2s 2 w . (5.13) 

We find (q = d, s) then 

U cS (b -> sf£) = U cS (b sj) - ^L^V*V tb £ [C^Qtin) + CKriQKfi)] (5.14) 

71 i=9,10 

where 

Q 9 = (s lfl P L b)(£Y£), Qio = (s llA P L b)(£Y^), (5.15) 

Q' 9 = (s lfl P R b)(£^£), Q' 10 = {s llx P R b){£^ lh £) ■ (5.16) 
For the Wilson coefficients we find 

1 1 A sb (Z) A^(Z) 

sin 2 6 W C 9 = [ VY Y (x t )-4sm 2 6 w Z (x t )} - — — j L \/ v 1 \ (5.17) 

S'SM M Z V ts V tb 

sin 2 9 W C W = - VY Y ( Xt ) - ^^ A %ly T(Z \ (5.18) 

fi'SM M Z V ts V tb 

sin 2 9 C- 11 A *W A vW ( 51 o) 

8 wCg ~ S&M v^v tb ' (5 ' 19) 

sin 2 C ~ 1 1 A ^ A ^) f5 20) 

sin Vw^io ~ --q-JJT \7*v ' ( 5 -^ 

9SM M Z V ts V tb 



where we have defined 

Af(Z)=Af(Z) + A>f(Z), 
Af(Z) = Af(Z)-Af{Z). 



(5-21) 



Here Y (x t ) and Z (x t ) are one-loop functions, analogous to X (x t ), that result from 
various penguin and box diagrams and given as follows 

^ r , \ Xt ( x t — 4 3x t log x t \ i \ 

^) = 7fc + (^J (5 - 22) 
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Figure 4: Tree level contributions of Z to the b — > s£ + £ effective Hamiltonian. 

„ , s 1 , 18xf - 163x, 3 + 259x, 2 - 108^ 32x? - 38xf - 15s 2 + \Sx t , 

Z»(*<) = - ? l°g*. + 14 4(x,-l)» 72(1,-1)* l0gI ' ' 

(5.23) 

r/y = 1.012 is QCD correction evaluated for m t = m t (m t ) [38| ll3]. 
Defining 

(5 - 24) 



Ag(Z)Aff(Z) 

1 2A%*(Z)Aff(Z ) 

4sin 2 ^ i^Mpp^ 



Z « = ^o(^) + - - J \' L " ' , (5.26) 



_ 1 2Aff(Z)Ag(Z) 

9 4sin 2 ^ &| M Mly t *y t6 



(5.27) 



we can write the Wilson coefficients as 

sin 2 9 W C 9 = Y q -A sin 2 6 w Z q , (5.28) 
sin 2 6 w C' g = Y' q - 4 sin 2 w Zj , (5.29) 
sin 2 fl^Cio = -K, , (5.30) 
sin 2 9 W C[ = -Y' q . (5.31) 

(5.32) 

The effective Hamiltonian for s — > d£ + £~~ transition can be obtained directly from [37] or 
from formulae given above by replacing q by K and appropriately changing the flavour 
indices. 
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6 Rare Decays 



6.1 K + — > 7r + vv and K L — > 7r°vv> 

Having at hand the effective Hamiltonian for s — > dvv transitions derived in Section 
it is now straightforward to obtain explicit expressions for the branching ratios B(K + — > 
7i + uu) and B(Kl — > tx q vv). Reviews of these two decays can be found in [HHI6]. 

The branching ratios for the two K — > ttuu modes that follow from the Hamiltonian in 
Section 15.21 can be written generally as 



B{K + -> 7i + uu) 



( ImX eff 

V ^ 5 



+ 



( ReX cS 



Pc{X) 



B{K L 7T°Z/Z/) = K L 



( ImX eff 



where [17] 



K+ = (5.36 ± 0.026) ■ 10" 11 , k l = (2.31 ± 0.01) ■ 10 



-10 



and [39H3I1I1ZISH]- 

P C {X) = 0.42 ± 0.03. 
The short distance contributions are described by 

X cS = V*V U (X L (K) + X R (K)), 



(6.1) 
(6.2) 

(6.3) 
(6.4) 

(6.5) 



where Xl,r(K) are given in (15. 7h and (15. 8p . They include both the SM contributions 
from Z-penguin and box diagrams and the tree-level Z contributions calculated here. 
The numerical analysis of both decays is presented in Section |HJ 

Experimentally we have jl9] 



+U.5N . in- 
10.5J iu 



B(K + -> ti + uu) cxp = (17.3 
and the 90% C.L. upper bound [50] 

B(K L 7r°z/z/) cxp < 2.6 ■ 10" 8 



ii 



(6.6) 



(6.7) 



In the SM one finds HUHZ 



SM 



(8.5 ±0.7) • 10 



-ii 



B(K L 7t°z/z7)sm = (2.6 ± 0.4) • 10 



-ii 



(6.8) 
(6.9) 
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where the errors are dominated by CKM uncertainties. This should be compared with 
the experimental values given in ( 16. 6 j) and (I6.7p . Clearly we have to wait for improved 
data. 

It should be emphasized that B{Kl —> k vv)sm depends sensitively on CKM parameters, 
in particular not only on \V c b\, as is also the case of B(K + — > n + uu)sM, but also on \V u b\ 
as it is a CP-violating observable. In our case the optimal value of \V u b\ will be 0.0037 
implying the central value for this branching ratio at 2.9- 10~ n , consistent with the value 
quoted above but higher. On the other hand our central value for B(K + — > 7T + uu)su 
agrees well with the one given in [12] . 

6.1.1 K, tt°£+£- 

The rare decays Kl — > n°e + e~ and Kl — > 7r°/i + / u _ are dominated by CP-violating 
contributions. The indirect CP-violating contributions are determined by the measured 
decays K$ — > tt°£ + £~ and the parameter Ek in a model independent manner. It is the 
dominant contribution within the SM where one finds [51] 

nn 11 , (6.io) 
nn 11 , (6.ii) 

with the values in parentheses corresponding to the destructive interference between 
directly and indirectly CP-violating contributions. The last discussion of the theoretical 
status of this interference sign can be found in [52] where the results of [53TI55] are 
critically analysed. From this discussion, constructive interference seems to be favoured 
though more work is necessary. In view of significant uncertainties in the SM prediction 
we will mostly use these decays to test whether the correlations of them with Kl — > iftvv 
and K + — > 7i + uu decays can have an impact on the latter. To this end we will confine our 
analysis to the case of the constructive interference between the directly and indirectly 
CP-violating contributions. 

The present experimental bounds 

B(K L -> 7r°e + e~) oxp < 28 • 10" 11 [56] , B[K L -)■ 7rVV~)ex P < 38 • 10~ n [57] , 

(6.12) 

are still by one order of magnitude larger than the SM predictions, leaving thereby large 
room for NP contributions. In the LHT model the branching ratios for both decays can 
be enhanced at most by a factor of 1.5 j58j|59]. Slightly larger effects are still allowed 
in RSc [37]. Much larger effects have been found in general Z' models [3]. However our 



B(K L -> ttVOsm = 3.5418* (1.5618S) • 
B(K L AV)w = 1-41^:1 (0.95i°;g) • 
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numerical analysis demonstrates that NP effects in these decays in the TUM are even 
smaller than in LHT. 



In the LHT model, where only SM operators are present the effects of NP can be com- 
pactly summarized by generalization of the real SM functions Yo(xt) and Zo(xt) to two 
complex functions Yk and Zk, respectively. As demonstrated in the context of the cor- 
responding analysis within RSc [37], also in the presence of RH currents two complex 
functions Yk and Zk are sufficient to describe jointly the SM and NP contributions. 
Consequently the LHT formulae (8.1)-(8.8) of [58] with Yk and Zk given below can be 
used to study these decays in the context of tree-level Z exchanges. Application of these 
formulae for general Z and Z 1 exchanges can be found in [3j. The original papers behind 
these formulae can be found in [51 ][53l, [5^ | [60 ] [6T] . 

Using the formulae of [37J we find 



Yi 



K 



r) Y Y (x t ) + 



Z K = Z (x t ) + 



1 



Af(g) 



A s v d {Z) 



V t * s V td 



4 sin 2 6 W 



2Af(Z) 
. M 2 zl g 2 SM 



A sd 



(Z) 



v t * s v u 



(6.13) 
(6.14) 



where Ay is defined in (15.211) . 



6.2 K L n + V 



As discussed in [37] in models with tree-level gauge boson exchanges in this decay the 
real function Y (x t ) is replaced by the complex function 



Y A (K) = r, Y Y (x t ) + 



[Af(Z)] 



Af(Z) - Ag(g) 
V t *Vta 



\Y A (K)\e ie r. 



(6.15) 



Only the so-called short distance (SD) part to a dispersive contribution to Kl —> n + n 
can be reliably calculated. We have then following [62] (A = 0.2252) 



in 2 



B(K L fi + fi~) S D = 2.03 • 10- 9 [P c (Y K ) + A 2 R t \Y A (K)\ cos^] 



(6.16) 



where R t and A are defined through \V t d\ = \V us \\V c b\Rt and \V c b\ = AX 2 , respectively. 
Moreover 



j3* = 13-/3,-6$, P C (Y K ) = [!-'—) P C (Y K ) , (6.17) 



l-y)^(^) , 

with P c (Yk) = 0.113 ± 0.017 [63]. Here (3 and (3 S are the phases of V t d and V ts defined 
in (fODjl . 
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The extraction of the short distance part from the data is subject to considerable uncer- 
tainties. The most recent estimate gives [61] 

B{K L -> n V)sd < 2.5 ■ 10- 9 , (6.18) 

to be compared with (0.8 ±0.1) • 10 -9 in the SM [63]. The numerical results are discussed 
in Section[8l In fact we will find that this decay plays a very significant role in our analysis 
as was already signaled by the rough estimates presented in pQ. 



6.3 B 



d,s 



We will next consider the important decays B^ s — > that suffer from helicity sup- 

pression in the SM. This suppression cannot be removed through the tree level exchange 
of Z boson but in principle could be removed through tree level exchanges of the Higgs 
boson. However the flavour conserving Hfift vertex is proportional to the muon mass 
and in contrast to SUSY and general two Higgs doublet models this suppression can- 
not be canceled by a large tan enhancement. Therefore in what follows we restrict 
our attention to the contributions of the Z° boson both through SM penguin and box 
contributions and its generated tree-level exchanges, calculated in Section 15.41 

Following (37] and assuming that the CKM parameters have been determined indepen- 
dently of NP and are universal we find 



B(Bg + \l 

B(B q -> /i+yU" 



where Yj^{B q ) is given by 
Y A (B q ) = r] Y Y (x t ) + 



At(Z) - Ag(Z) 



Ya(B„ 



r] Y Y (x t ) 



(6.19) 



qb, 



v t *v tb 



\Ya{B q 



\e< q 



(6.20) 



As stressed in [65T467] §, when comparing the theoretical branching ratio B(B S — > yU + /i _ ) 
with experimental data quoted by LHCb, ATLAS and CMS, a correction factor has to be 
included which takes care of Ar s effects that influence the extraction of this branching 
ratio from the data: 



Here 



B(B S ) th = r(y s ) B(B S ) cxp , r(0) = 1. 

i-y 2 s 



r(y s 



1 + J^^pUs 



1 - A A AT y s 



(6.21) 
(6.22) 



5 We follow here presentation and notations of [66 ( 67 
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with 

Ar 

y s = t Bs = 0.088 ± 0.014. (6.23) 
The quantity *44r is discussed below. 

The branching ratios B{B q — >• are only sensitive to the absolute value of Ya(B q ). 

However, as pointed out in [6"6"|l67] in the flavour precision era these decays could allow to 
get also some information on the phase of Y^Bq). The authors of |67tl68] provide general 
expressions for A\ T and S s + as functions of Wilson coefficients involved. Using these 
formulae we find in our model very simple formulae that reflect the fact that Z° and not 
scalar operators dominate NP contributions: 

A X AT = cos(20?' - 2 VBs ), S; +fl . = sin(20*' - 2ip B ,). (6.24) 

Both A\ T and S s + are theoretically clean observables. 

While Ar^ is very small and yd can be set to zero, in the case of Bd — > one can 

still consider the CP asymmetry [HE], for which we simply find 

S d ^ =sin(2^-2^J. (6.25) 

As demonstrated in Section [8] the TUM model considered by us turns out to give values 
for B(B s d — > that uniquely differ from SM predictions. However, similarly to the 

SM and CMFV models we find in TUM to a very good accuracy 

A x Ar = 1, S'^. = 0, r{y s ) = 0.912 ± 0.014 (6.26) 

basically independent of NP parameters considered. These are definite predictions of 
TUM which will be tested one day. 

In our numerical analysis in Section [8] we will discuss then only the branching ratios for 
these decays. Now, the most recent results from LHCb read [21169] 

B(B S -»■ fj,+ fT) = (3.2±J:|) x HT 9 , B(B S //+//~) SM = (3.23 ± 0.27) x 10~ 9 , (6.27) 

B(B d ->■ fi + (T) < 9.4 x 10~ 10 , B(B d SM = (1-07 ±0.10) x 10" 10 . (6.28) 

We have shown also SM predictions for these observables [70] that do not include the 
correction r(y s ). If this factor is included one finds [661167] 

B(B S ->■ = (3-5 ± 0.3) • 10~ 9 . (6.29) 

It is this branching that should be compared in such a case with the results of LHCb 
given above. For the latest discussions of these issues see [66ti68ll70] . 
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As in TUM A\ r = 1 independently of NP parameters we will include the correction in 
question in the experimental branching ratio using the values in (16. 26 p . If this is done 
the experimental results in (16. 2 7ft is reduced by 9% and we find 

B{B S -> /i + /Ocarr = (2.9±i;f) x 10" 9 , (6.30) 

that should be compared with the SM result in (16.271) . While the central theoretical 
value agrees very well with experiment, the large experimental error still allows for NP 
contributions. In our plots we will show the result in (16.3QH . 

This completes the analytic analysis of the B s ^ —> ^ + ^~ decays. The numerical results 
are discussed in Section El 



6.4 B -> {X s ,K,K*}vv> 

Following the analysis of [71], the branching ratios of the B — > {X s , K, K*}uu modes in 
the presence of RH currents can be written as follows 

B{B -> Kvv) = B(B -> Kvv) m X [1 - 2r?] e 2 , (6.31) 

B(B K*vv) = B(B -»■ K*vv) su x [1 + I.3I77] e 2 , (6.32) 

B(B X a vv) = B[B -»■ X s z/z/) SM x [1 + 0.09r/] e 2 , (6.33) 

where 

2 _ \X L (B s )f + \X R (B S )\* ^ = -Re(X L (B s )X^(B s )) ^ 



\VxX (xt)\* ' ' \X L (B S W + \X R (B S )\* 
with XL n(B a ) defined in (15. lip . The issue of long-distance contributions to these short- 
distance formulae is discussed in [72] . 

The predictions for the SM branching ratios are [7TH73] 

B(B -»■ Kvv) SM = (3.64 ± 0.47) x 10" 6 , 
B(B -> KVp) sm = (7.2 ± 1.1) x 10" 6 , 

B{B -»■ X s z/z/) SM = (2.7 ± 0.2) x 10" 5 , (6.35) 
are respectively by factors of four, eleven and twenty below the experimental bounds [74T 



We would like already announce at this point that the TUM model satisfies easily the 
present experimental bounds and makes definite predictions for e and rj defined in (I6.34p : 

e>l, 77 «0, (6.36) 

The last result is the consequence of the strong suppression of right-handed contributions 
when all constraints are taken into account. More details will be given in Section El 
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6.5 B + -> t+v 



6.5.1 Standard Model Results 

We now look at the tree-level decay B + — > t + v which in the SM is mediated by the W ± 
exchange with the resulting branching ratio given by 

B(B+ -+ r+i/)sM = Glm ^ ml (l - ^f) 2 Fl + \V ub \*r B+ . (6.37) 

87T y m B+ J 

Evidently this result is subject to significant parametric uncertainties induced in ( 16. 37ft 
by Fb+ and V u b- However, it is expected that these uncertainties will be eliminated in 
this decade and a precise prediction will be possible. Anticipating this we will present 
the results for fixed values of these parameters. 

In the literature in order to find the SM prediction for this branching ratio one eliminates 
these uncertainties by using AMj, AMd/ AM S and S^k s [771175] and taking experimental 
values for these three quantities. This strategy has a weak point as the experimental 
values of AM djS used in this strategy may not be the one corresponding to the true value 
of the SM. However, proceeding in this manner one finds [78] 

B(B + -»■ r + z/) SM = (0.80 ± 0.12) x 10~ 4 , (6.38) 

with a similar result obtained by the UTfit collaboration [77] and CKM-fitters. 

For quite some time this result was by a factor of two below the data from Belle and 
BaBar. However this disagreement of the data with the SM softened significantly with 
the new data from Belle Collaboration [16]. The new world average provided by the 
UTfit collaboration [17] 

B(B + -> r+z/) exp = (0.99 ± 0.25) x 10~ 4 (6.39) 

is in perfect agreement with the SM, even if NP providing a slight enhancement of this 
branching ratio is presently favoured. 

The full clarification of the left room for NP in this decay will have to wait for the data 
from Super-B machine at KEK. In the meantime hopefully improved values for Fg+ 
from lattice and \V u j)\ from tree level decays will allow us to make a precise prediction for 
this decay without using the experimental value for AMd- In TUM it will turn out that 
the favoured value of \V u b\ = 0.0037 in this model implies through (16. 37ft central value 
for B{B + — > t + v) pd 0.88 x 10 -4 , that is very close to the data and the question arises 
whether modification of W ± couplings in TUM still allows to keep this agreement. 
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6.5.2 Effect of Modified W ± Couplings 



In the presence of modified W couplings we find 



B(B + -> t + v) 
with 



64vrM^ 



m B +m T 1 



m 



2 \ 2 



m, 



Fl + T B+ \AY{^t-^ U R)\\ (6-40) 



Al b = -j=V ub + Al 3 (W), Ag = Ag(W0, A " = -7f 



(6.41) 



with A] : 3 R (W / ) given in Section E2] For Af = Af = this expression reduces to 
the SM expression in ( 16.37p . Using our estimates of the corrections to W couplings in 
subsection 13.101 we find 



V2 



'(tf) 



< 10 



-2 



(6.42) 



and consequently NP corrections to B(B + — >■ r + z/) can be safely neglected. 



7 Basic Structure of New Physics Contributions in 
the TUM 

7.1 Preliminaries 

The spirit of our analysis presented in this paper differs significantly from the one of our 
first paper pQ and many papers found in the literature in which the main goal is to find 
out whether a given model is roughly consistent with the present experimental bounds. 
In doing this quite often one imposes the constraint that NP contributions are at most 
as large as the SM contributions. 

In view of increased precision in experimental data and in the theory to be expected 
in this decade such a passive approach to the tests of new model constructions is not 
satisfactory. Below we will follow a more aggressive approach by identifying correlations 
between various observables predicted by the TUM and proposing various tests of this 
scenario. 

In many NP extensions of the SM, like LHT model, RS scenarios of various sort and 
supersymmetric models the identification of specific correlations between various ob- 
servables is challenged by the multitude of new free parameters present in these models. 
However, in the case of TUM we are in a comfortable situation that our model does 
not have many free parameters to describe FCNC processes as several of its fundamental 
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parameters have been already used to describe successfully the spectrum of quark masses 
and the values of the CKM parameters. In fact as we already discussed previously once 
the common mass of the vectorial heavy fermions has been fixed, only three real and 
positive definite parameters are left to our disposal and NP contributions to all FCNC 
processes in the down quark sector are entirely given in terms of these parameters and 
the parameters of the SM. Very important is the fact that the sole CP-violating phase 
in the TUM equals the KM phase and is equal to the angle 7 in the unitarity triangle. 

7.2 Facing the Anomalies in AF = 2 Data 

Before entering the details let us first ask the question how TUM faces the anomalies in 
the data for AF = 2 observables identified first in [32],[T9] . Indeed the SM does not offer 
a fully satisfactory description of these data. Here the prominent role is played by the 
cr — S^Ks tension within the SM. In this context it should be emphasized that because of 
this tension the pattern of deviations from SM expectations for other observables depends 
often on whether ex or S^k s is used as a basic observable to fit the CKM parameters. As 
both observables can in principle receive important contributions from NP, none of them 
is optimal for this goal. The solution to this problem will be solved one day by measuring 
the CKM parameters with the help of tree-level decays. Unfortunately, the tension 
between the inclusive and exclusive determinations of \V u b\ and the poor knowledge of 
the angle 7 from tree-level decays preclude this solution at present. However a good 
agreement of the SM value for the ratio AM S / AM^ with the data and the fact that the 
present determinations of 7 are in the ballpark of 70°, imply two basic scenarios for the 
three observables |T4&|, £r and S^k s - Moreover for fixed 7 and AM S / AM d the latter 
three observables are strongly correlated within the SM with each other. We have: 

• Exclusive (small) \V U \,\ Scenario 1: \sk\ is smaller than its experimental deter- 
mination, while S^Ks is very close to the central experimental value. 

• Inclusive (large) \V u b\ Scenario 2: \ek\ is consistent with its experimental de- 
termination, while S^k s is significantly higher than its experimental value. 

Thus dependently which scenario is considered we need either constructive NP contribu- 
tions to \ek\ (Scenario 1) or destructive NP contributions to S^k s (Scenario 2). However 
this NP should not spoil the agreement with the data for S^k s (Scenario 1) and for \ek\ 
(Scenario 2). 

While introducing these two scenarios, one should emphasize the following difference 
between them. In Scenario 1, the central value of \ek\ is visibly smaller than the very 
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precise data but the still significant parametric uncertainty due to \V c b\ 4 dependence in 
\ek\ and a large uncertainty in the charm contribution found at the NNLO level in [21] 
does not make this problem as pronounced as this is the case of Scenario 2, where large 
| V u b | implies definitely a value of S^k s that is by 3a above the data. 

Now models with many new parameters can face successfully both scenarios removing 
the deviations from the data for certain range of their parameters but in simpler models 
often only one scenario can be admitted. For instance in models with constrained MFV 
there are no NP contributions to S^k s but NP contributions to Ek are possible so that 
Scenario 1 is selected in that framework. On the other hand in 2HDMj^py [80] there 
are no relevant NP contributions to Ek but the presence of flavour blind phases can 
have a significant impact on S^k s so that Scenario 2 is selected. On the other hand in 
models with U(2) 3 flavour symmetry both scenarios can be accommodated implying an 
interesting triple correlation between the values of S^k s , and \V u b\ [81 J. 

As we will demonstrate in the next subsection first at a semi- quantitative level, the TUM 
cannot improve the description of AF = 2 data relative to the SM for the following 
reasons: 

• NP contributions to B° s d — B° sd mixings are very small so that 

S^Kg, Sm, AM d , AM S (7.1) 

remain SM-like, although the predicted value of S^k s depends sensitively on the 
chosen value of \V u b\- We simply find that within 1% accuracy 

C Bd = C Bs = 1, ip Bd = VBs = 0. (7.2) 

• NP effects in Ek can in principle be large so that at first sight one could choose 
exclusive scenario for |V^&| and enhance Ek through NP contributions. However in 
TUM the LR contributions uniquely suppress Ek relatively to its SM value and have 
to be compensated by LL contributions. But the increase of these contributions is 
bounded from above by the upper bound on Kl —> and we find that in order 
not to violate this bound LL contributions can just compensate LR contributions 
so that Ek is very close to its SM value. 

In view of this situation, within the TUM the optimal choice for \V u b\ is between the 
Scenarios 1 and 2. Inspecting then the dependence on 7 we are then lead to the following 
optimal central values for the four CKM parameters that we will use in our analysis: 

\V US \ = 0.2252, \Vcb\ = 0.0406, \V ub \ = 0.0037, 7 = 68°, (7.3) 
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TUM 


Experiment 


\ £ k\ 


2.00(22) • 10~ 3 


2.228(11) x 10" 3 


SipK s 


0.725(25) 


0.679(20) 




0.039 


0.002(9) 


AM S [ps- 1 ] 


19.0(21) 


17.73(5) 


AM, [ps- 1 ] 


0.55(6) 


0.507(4) 


B(B + -> r+z/ T ) 


0.88(14) • 10" 4 


0.99(25) x 10" 4 



Table 3: TUM predictions for various observables for \V u b\ = 3.7 • 10 3 and 7 = 68° 
compared to experiment. 



where the values for \V US \ and \V c b\ have been measured in tree level decays. Moreover 
the value for 7 is consistent with the ratio AM^/AM S in the model considered: 

(AMA = (AM,\ =345±30 

In Table [3] we summarize for completeness the results for various observables in the TUM 
model that for the input in (17.31) and in Table HI are the same as in the SM. 

We note that: 

• \ek\ is lower than the central experimental value but still within la when hadronic 
and parametric uncertainties are taken into account. 

• SWs * s by 2c larger than its central experimental value. 

• is very close to the experimental value. 

• AM, and AM,, although slightly above the data, are both in good agreement with 
the latter when hadronic uncertainties are taken into account. In particular their 
ratio is in an excellent agreement with data. 

• Choosing higher value of 7 would bring \e K \ closer to the data. But then also AM, 
would be larger implying worse agreement with the data for AM, and AM S /AM,. 

These results depend on the lattice input and in the case of AM, on the value of 7. 
Therefore to get a better insight both lattice input and the tree level determination of 7 
have to improve. 
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Thus the description of AF = 2 observables in TUM, similarly to the SM, is not fully 
satisfactory but sufficiently good that we can continue our analysis and investigate the 
predictions in this model for 

B Sjd ^fi + fi~, B — >• K*(K)£ + £~ , B->X a £r (7.5) 

B^X s uu, ,B^K*vv, B^Kvv, (7.6) 
K+^n+uu, K L ^7T°uu, K L -+n + fj,-, K L ^7T°e + £-. (7.7) 

Now, it is known from [3] that in the case of tree-level FCNCs mediated by Z° generally 
NP effects in AF = 1 processes are larger than in AF = 2 transitions, while the opposite 
is true for Z' tree-level exchanges with Mz> of order few TeV. In the concrete model 
considered here there is an additional suppression factor that makes NP effects in AF = 2 
transitions to be much smaller than in AF = 1 transitions. As these properties are very 
important for our analysis we explain them here in explicit terms: 

• Considering first as in [3] a general tree-level neutral gauge boson (A) contributions 
to FCNC processes, a tree-level A contribution to AF = 2 observables depends 
quadratically on A l [ R (A) / Ma, where A l [ R (A) are flavour-violating couplings with 
i, j denoting quark flavours. For any high value of Ma, even beyond the reach of the 
LHC, it is possible to find couplings A l [ R (A) which are not only consistent with the 
existing data but can even remove certain tensions found within the SM in AF = 2 
processes. The larger Ma, the larger A % [ R (A) are allowed: A % [ R (A) ~ a^M^ with 
ciij adjusted to agree with AF = 2 data. Once A l [ R (A) are fixed in this manner, 
they can be used to predict A effects in AF = 1 observables. However here 
NP contributions to the amplitudes are proportional to A l [ R {A)/M\ and with the 
couplings proportional to Ma, A contributions to AF = 1 observables increase with 
decreasing Ma without changing AF = 2 transitions. Eventually, for sufficiently 
low values of Ma the bounds on AF = 1 processes become stronger than the 
constraints on AF = 2 processes requiring the coefficients to be smaller than 
obtained from AF = 2 constraints alone. Therefore in turn for these low masses of 
Ma, as is the case of Mz, even if NP effects mediated by Z° in AF = 2 transitions 
are small, significant effects in rare K and B decays can be found. In this manner 
flavour-violating Z couplings turn out to be an important portal to NP, in our case 
the physics of vectorial fermions. 

• In the specific model considered by us this disparity between AF = 2 and AF = 1 
transitions is enhanced by the fact that the coefficients in the couplings A l [ R (Z) 
are proportional to v 2 /M 2 and consequently strongly suppressed. As AF = 2 
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amplitudes are quadratic in these couplings and the interference between SM and 
NP contributions dominating NP effects in rare decay branching ratios is linear in 
them, the disparity in question is in our model much larger than could be expected 
on the basis of the general analysis in [3]. 

Now, even if NP contributions to B® d — B Q sd mixings are found to be very small, the 
structure of NP contributions to Ek forces flavour violating A| d (Z) couplings to be 
sufficiently large in order to compensate the negative contributions of LR operators. 
But as A S ^(Z) couplings are bounded from above by K L — > this compensation 

is only possible if RH couplings are suppressed. Recall that hadronic matrix elements 
of LR operators are in the K system not only enhanced through renormalization group 
effects but also receive a large chiral enhancement. 

Thus at the end NP contributions to rare K decays are governed by LH couplings of Z° 
to quarks with RH couplings being subleading. Due to small number of free parameters 
in TUM this structure is transferred to rare B decays. At the end the predictions of 
TUM for rare K and B decays are forced to differ from SM predictions and the deviations 
from the SM take place in a correlated manner. In order to have a better insight in this 
structure we will now derive approximate expressions for various amplitudes and master 
functions governing various observables so that the pattern just discussed will be seen 
in explicit terms. 



7.3 Structure of Z Q Couplings 

The starting point of our discussion are suitable approximations for the A®, which are 
sums of three different terms (summation over k = 1,2,3) 

~ a d ._. v 2 Q2 ~ D _ mi m 3 v k \v k] 



We find 



A?^^(e?X t+ e?^) A^^(^ + i ^) (7.9) 



v c Q2 x d ad „ m bm d ( X u ReAg 

M 2 £ 3 At A R31 ~ ^ ^ 2 + ^ 



A D ^^-e Q2 \ s A D ~ msmh ( K I ReK ) (7 11) 

A L32~ n/r o e 3 At A R32~ nj2 \ Q2 + Q2 ' 



M 2 1 m2 M 2 \,? 2 £ 2 

where \ and are defined in H4.3[) 
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Already at this stage we observe that the LH couplings in B% and B® systems are of 
MFV type, while in the K system this property is broken by the term involving A c . The 
RH couplings are suppressed through quark masses but in order to estimate their size 
Xij have to be constrained through FCNC processes as done below. 

Using then ( 13. 5 7ft we can express these couplings entirely in terms of €® , Xi 3 and X 2 3, 
quark masses and the elements of the CKM matrix. We find 

A?, « (A. + W ReA c ) i», w ™ £ + ^) 

(7.12) 

~ a d „ v \Q2 x d ad „ m bm d 1 / Xj Re\ d c 

A L31 ~ T72 e 3 A t AR31 r ~ ™ I ~ ,o T 



M 2 3 « ™ M 2 £ f V*ulK»| 2 A 2 3 |U cfe | 2 



(7.13) 



v2 m,s W „ m sm b 1 / A* ReA;! 

M^ 3 A * Ar32 ~ M 2 £ f Ui 2 3 |K 6 | 2 + X^I^P 



(7.14) 

We recall that the masses of quarks entering these expressions are evaluated at the high 
scale M and for M = 3 TeV take the values 

m d (M) = 2.3 MeV, m s (M) = 45 MeV, m b (M) = 2.4 GeV . (7.15) 



7.4 Structure of AF = 2 Amplitudes 



Using these results we find the tree-level Z contributions to mixing amplitudes Mf^ and 
Mf 2 . We list them in Appendix O Inserting in the latter formulae our nominal values 
of the CKM parameters in (17.31) and quark masses in (I7.15P we obtain the following 
approximate results for quantities of direct interest: 



Im(Mg) z 
Im(M*) SM 

Re(M^) z 
Re {M K ' 



i2;sm 



4.0 



3.0 



/3TeVy 
V M ) 

/3TeV\ 4 



e? 4 (l + 1.2 X. 



2 1 

2.3 ) 



0.18 



\P™{K)\ 



-Q4 y-4 

'3 ^23 



0.30^^(^)1 



X 13 



Xno 1 

-f + 0.86- 

A 13 



X13 



(7.16) 

10" 5 
(7.17) 



Re 



M 



12 



[m- 



2.6 



/3TeVY 



-Q4 



10" 



(7.18) 
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Im 



M 



SM 



Re- 



Im 



-3.0 



2.5 



-^13 



2.3 



/3TeV\ 4 

U j 

/3TeV\ 4 



-10" 



-Q4 



4.2 



^13 



-10" 



10" 



« 2 S )SM 

For M = 3 TeV we find 
\P™{K)\ 



\ M 



39, 



\P£ R (B a 



-10" 



(7.19) 

(7.20) 
(7.21) 



4.4, 



P! LR (5 S 



4.5. 



(7.22) 



All P 1 " 11 are negative and the minus signs have been included in the formulae above. 
As seen in Appendix O the above expressions can allow us to estimate the size of new 
contributions to ex, SAMk, Cb„, SW= and S, 



The last four results signal that NP effects in B® — B® systems are suppressed. However, 
in order to prove it we have to know the size of X 13 . £3 is 0(1). To this end it is useful 
to look first at AF = 1 transitions. 



7.5 Structure of AF = 1 Amplitudes 

Proceeding in the same manner we obtain approximate expressions for the relevant 
quantities in AF = 1 observables which we list in Appendix [D] Inserting the nominal 
values of CKM parameters and quark masses we find 

AX L (K) = X L (K) - X(x t ) = 0.31 f^-Yef [l + l.lX&eW-M] (7.23) 

/3TeV\ 2 e^-M 

AX L (B d ) = 0.3l(^)V (7.25) 
X^)^-,7(^) 2 -i^e-^)l0-3 (7.26) 

AX L {B S ) = 0.31 ( ^ ) ef (7.27) 



X R (B S ) = 0.31 



3TeV x 




M , 


3TeV x 




M , 


) ~k 
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5.5e" i7 2.2 



+ 



Y 2 Y 2 
^13 ^23 



10~ 3 (7.28) 



Moreover, we find 

AY A {K) = AX L (K) - X R (K), AY d = AY S = AX L (B d ) = AX L (B S ) (7.29) 
with other relations listed in Appendix [D] AY q are NP corrections to Y q in (I5.24p . 

7.6 The Interplay of AF = 2 and AF = 1 Transitions 

With these formulae at hand we can now understand the pattern of NP effects that we 
have outlined in Subsection 17.21 Indeed: 

• As seen in ( 17. 161) the last term representing LR contribution to ek and being en- 
hanced through |P! LR (^)| « 39 suppresses \ek\ instead of enhancing it as required 
by the data. For M < 3 TeV this contribution is a problem for TUM. 

• The solution to this problem is a sufficiently small value of 1/ Xi% accompanied by 
sufficiently large values of and X 2 3- However < 1. 

• The required suppression of I/X13 because of Ek and the small quark masses 
multiplying it in the expressions for corrections to the RH master functions as 
given in in Appendix [D] imply that AF = 1 transitions are dominated by LH 
currents. This is explicitly seen in (17.241) . ( 17.261) and ( 17.281) . 

• On the other hand e® and X23 are bounded from above by BiK^ — > / u + fi~) as with 
the suppression of RH currents AYa(K) AXl(K) and with increasing and 
X23 the upper bound in ( 16.181) is approached. However, as our numerical analysis 
in the next section shows, this bound still allows for values of and X23 necessary 
to compensate the LR contribution to Ek in ( I7.16p . But then the net effect of NP 
in Ek is very small and \ V u b\ has to be larger than its exclusive determinations in 
order for Ek to be within la from the data. 

• Remarkable are also results in (I7.23p . (I7.25P and (I7.27P which imply that all rare 
decay K and B q branching ratios considered by us are predicted to be enhanced 
over their SM values. Moreover, the enhancements of branching ratios for B s 
and Bd rare decays satisfy CMFV relations. On the other hand the usual CMFV 
relation between rare B q decays and rare K decays is broken by non-vanishing A23 
but this effect, as explained in the next section is only present in K + — > -k + vv. 

With this general view in mind we can now enter the numerical analysis. 
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G F = 1.16637(1) x l(T 5 GeV- 
M w = 80.385(15) GeV 
sm 2 6 w = 0.23116(13) 
a(M z ) = 1/127.9 
a s (M z ) = 0.1184(7) 



m u (2GeV) = (2.1 ±0.1)MeV 
m d {2 GeV) = (4.73 ± 0.12) MeV 

m s (2 GeV) = (93.4 ± 1.1) MeV 

m c {m c ) = (1.279 ± 0.013) GeV 
m b (m b ) = 4.19±g;S| GeV 
m t (m t ) = 163(1) GeV 
M t = 172.9 ±0.6 ±0.9 GeV 



m K = 497.614(24) MeV 

F K = 156.1(11) MeV 

B K = 0.767(10) 

n e = 0.94(2) 

771 = 1.87(76) 

m = 0.5765(65) 

773 = 0.496(47) 

AM K = 0.5292(9) x lO^ps -1 
lejcl = 2.228(11) x lO" 3 



[82J 
[82] 
1821 



[83J 
1831 



[82J 
E31 



[23] 
[H2] 

E21 



B(B -> X sl ) 
B(B+ -> r+i/) 



(3.55 ±0.24 ±0.09) x 10 
= (1.64 ±0.34) x 10- 4 



4 [82] 
E2] 



r B ± = (1641 ± 8) x IO- 3 



ps 



m Bd = 5279.5(3) MeV 
m Bs = 5366.3(6) MeV 
F Bd = (190.6 ± 4.6) MeV 
F Bs = (227.7 ± 6.2) MeV 
B Bd = 1.26(11) 
B Bs = 1.33(6) 
B B jB Bd = 1.05(7) 



F Bd JB Bd = 226(13) MeV 



F Bs \JB Bs = 279(13) MeV 
f = 1.237(32) 
r/B = 0.55(1) 
AM d = 0.507(4) ps" 1 
A M s = 17.73(5) ps" 1 
S fKs = 0.679(20) 
= 0.0002 ± 0.087 

t(B s ) = 1.471(25) ps 
r(B d ) = 1.518(7) ps 



L82J 

m 

[83] 
[83] 
E3 



[83] 

[83] 
[83] 

22] 

82] 
87] 
82] 



[89J 
1891 



|K* I = 0.2252(9) [82] 
\V cb \ = (40.6 ±1.3) x 10" 3 [82] 
= (4.27 ±0.38) x 10" 3 [82] 
|V; e 6 xcl | = (3.12 ±0.26) x 10" 3 [83] 



Table 4: Values of the experimental and theoretical quantities used as input parameters. 

8 Numerical Analysis 



8.1 Procedure in the Trivially Unitary Model 

It is not the goal of this section to present a full-fledged numerical analysis of the TUM 
including present theoretical and experimental uncertainties as this would only wash out 
the effects we want to emphasize. Therefore, in our numerical analysis we will choose 
the values in (17. 3p as nominal values for four CKM parameters. The remaining input 
parameters are collected in Table HJ In any case as NP effects in B^ s — B® s are very 
small, the hadronic uncertainties in this sector are identical to the ones in the SM. They 
are more important in the case of Ek, as they play the role in the compensation of LR 
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contributions by LL ones and consequently have an impact on the allowed values of 
X 13 and X 2 3 in (I3.57p . Also the uncertainties in Fb are relevant for the predictions of 
B(B S ^ but lattice calculations made significant progress in the last years [83J. 

The observables analysed numerically in this section depend on thirteen real parameters 
and one complex phase that is equal to the angle 7 in the unitarity triangle. In the 
TUM, once the six quark masses and the CKM parameters have been determined as 
explained above, the model contains four positive definite real parameters 

M, si a ? 3 , 4, (8.1) 

with sf 3 and S23 smaller than unity and 0.80 < e 3 < 1.0 because of the value of the 
top quark mass. For fixed M eliminating the 10 parameters 2 , £12 3; £ i2 35 ^ s i2 i R 
favour of 3 CKM mixing angles, 1 phase and 6 masses, one can find the couplings 
A % 1 R for W ± , Z° and H° as functions of s^, sf 3 and £3 . In what follows we will set 
M = 3.0 TeV which on one hand is still in the reach of the LHC and on the other 
hand is sufficiently large so that the upper bound on B{Kl — > and also the LHCb 

bounds on b — > s£ + £~ transitions are satisfied. As shown in [TJ for such values of M the 
model is also consistent with electroweak constraints. We can then vary sf 3 , and e 3 . 
All observables considered by us depend now on only three real and positive definite free 
parameters implying various correlations that we are going to expose below. 



8.2 Phenomenology of A.F = 2 observables 

As we already discussed in Section [7] NP effects in B® sd — B® sd mixings are negligible. 
However they can be in principle large in Ek- Yet requiring that§| 

0.75 < — . < 1.25, 2.0 x HT 3 < \e K \ < 2.5 x 10~ 3 (8.2) 
(AMx)sm 

and imposing the bound on B(Kl — > /i + /i~)sD in (16. 18ft we find, as seen in Fig. |5] (left), 
that \ek\ is forced to be at the lower end of the range in (18.21) and very close to the SM 
expectation for the CKM parameters chosen by us. 

The colour coding in this plot and the following plots is as follows: 

• Green range is allowed by Ek through (18.21) and Ki — > ^ through ( 16.181) . 

• Yellow range is allowed by ek but not Kl —> 

• Purple range is forbidden by ek but allowed by K L — > fi + fi~ . 



3 The ranges chosen in (|8.2j) indicate theoretical and parametric uncertainties in both observables. 
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BR(K,.-»/i + /0 & \e K \ 

10 1 1 .0 1 




1.0 1.5 2.0 2.5 3.0 0.0 0.2 0.4 0.6 0.8 1.0 

\e K \ [10~ 3 ] 4 

Figure 5: Left: B(Kl —> versus \ek\ for M = 3 TeV and \V U b\ = 0.0037. Green 

points are compatible with both bounds for \e k \ (EUP and B(K L — > A6.18\) , yellow 

is only compatible with \ek\ and purple only with B(Kl — > /i + /i~). The red point corre- 
sponds to the SM central value. Right: Allowed region in the parameter space (sf 3 ,S2 3 ) 
for e 3 = 0.9 due to \ek\ and B{Kl — > fi + fi~) bounds (same colour coding as left). 
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Figure 6: X 23 versus I/X13. Colour coding as in Fig. O 



• Points that are forbidden by both ek and Kl — >■ H + are not shown. 

The corresponding regions in the space (sf 3 , s^) are shown in Fig. [5] (right) for e 3 = 0.9. 
We observe that in the allowed region > 0.5, and sf 3 < 0.6. 

The fact that NP effects in Ek are very small implies the following approximate relation 
between Xi 3 , X23 and e 3 

Ef (1 + 1.2 X 2 2 3 ) - 0.18 |P V1 K)I « (8.3) 



^13 
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In Fig. [6] we show the four regions in the space (X23, 1/X13). We observe that the allowed 
ranges represented by the green region are roughly: 

0.4 < X 23 < 2.0, 0.03 < -J- < 0.70. (8.4) 

-^13 

We also recall that for M = 3 TeV in order to reproduce the top-quark mass we have 
0.80 < £3 < 1. Consequently we conclude that the effects of right-handed couplings 
in AF = 1 decays are negligible and NP contributions to these decays are dominated 
by left-handed Z couplings to quarks. This is in contrast to rare decays in RS model 
with custodial protection [37] where NP effects in these decays were governed by flavour- 
violating RH couplings of Z° to quarks. 



8.3 B(B d , s fj,+fi-) and AM d /AM s 

In models with CMFV these observables are related through a theoretically very clean 
relation [90J that in the MTFM and generally in models with non-MFV sources and new 
operators gets modified as follows: 

B(B S -> //+//-) B Bd r(B s ) AM S 
B(B d ^n+n-) B Bs r{B d )AM d r ' 

where 

= Y A (B 3 ) 
Y A (B d ) 

with r = 1 in CMFV models but generally different from unity. 

However in TUM when all constraints are taken into account we find r = 1 and the 
relation (18. 5p is satisfied very well. Moreover CB ds = 1- Yet, the values of B(B q — >• 
differ significantly from SM prediction. This we show in Fig. [7J indicating the 
experimental la and 2a ranges for B(B S — > [i + The striking prediction of TUM are 
uniquely enhanced values of both branching ratios, moreover in the allowed green region 
these enhancements take place in a correlated manner: the slope of the straight line in 
Fig. [TJis given by the formula (18. 5 1) with r = 1. We find the ranges: 

4.2xl0~ 9 < B(B S ->-/i + /0 <6.0xl0" 9 , 1.3x 10~ 10 < B(B d ->■ <2.0xl0" 10 

(8.7) 

where we added parametric uncertainties not shown in the plot. 
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(8.6) 
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B(B 8 ->/iV) [1CT 9 ] 



Figure 7: B(B d — )■ ) versus B(B S — »■ ). Colour coding as in Fig. [3 Cray 

region: exp. 1 and 2a range of B(B S — >■ fi + fi~) (see Eq. K6. 30\) ). 



8.4 The K + — )- 7r + z/z/ and A"/ — ► 7r°i/i/ Decays 

In order to understand the pattern of NP contributions to these decays we calculate X c g 
in ( 16. 5 p in the TUM. Neglecting the small RH contributions we find 



-X"efT — — |Vt s ||Vtd| 

where 



+ 0.37 (*™y e ?J& 



X( Xt )=X(x t )+0.31 (^) 2 ^f 



1.8) 



5.9) 



The first term in f)8.8p describes a typical CMFV contribution with a modified basic 
function X that is uniquely larger than X in the SM. This increase is governed by 
e® 2 /M 2 . The second term does not carry any new phases and goes beyond CMFV. It 
contributes only to K + — > 7i + uu and modifies the usual CMFV correlation between the 
branching ratios for these two decays. 

In Fig. [8] we show the correlation between B(K + — > 7i + uu) and B(Kl —> 7r°i/z/) in 
TUM. The experimental la-range for B(K + — > tt + vv) [19] and the model-independent 
Grossman-Nir (GN) bound [91] are also shown. We observe that B(Kl — > 7T uu) can 
be as large as 4.4 • 10~ n , that is roughly by a factor of 1.5 larger than its SM value 
(3.0±0.6)-10~ n , while being still consistent with the measured value for B(K + — > -k + vu). 
The latter branching ratio can be enhanced by at most a factor of 2 but this is sufficient 
to reach the central experimental value |49j in (16.61) . 
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Figure 8: B(Kl — > ti uu) versus B(K + — > 7i + uu). The black line on the left is the GN 
bound. Gray region: experimental range of B(K + — > -k + vu) (see Eq. $6.60 ). Colour 
coding as in Fig. 



The plot has a shape which differs from the one encountered in LHT model or Z' models. 
The expression in ( I8.8p explains what is going on: 

• For a fixed value of £^ the branching ratio B(K^ — > 7r°z/z/) is fixed, while B{K + — > 
-n + vu) depending in addition on X 2 3 can take a significant range of values bounded 
by Kl — » and We have then a straight horizontal line. This line moves 
up with increasing £3 . 

• For fixed M both branching ratios increase with increasing e®. 

Thus if X 2 3 would vanish, we would just have the case of CMFV. Varying £3 we would 
get a straight line on which both branching ratios would increase with increasing £3 . 
This is the line which fully describes BiK^ — > 7r°z/p) in the TUM. However, X 2 3 cannot 
vanish and is bounded from below in order to balance negative contributions from LR 
operators to Ek- This implies the shape in the Fig. |HJ 

8.5 Correlation between Kl —± and K + —± tt + vv 

Next in Fig. [9]we show the correlation between BiK^ — > /U~V~)sd an d B(K + — > tt + vv). 
As both decays are CP-conserving, a non-trivial correlation is generally expected. The 
following observations should be made: 
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Figure 9: BiK^ — > versus B(K + — > tt + vv). Colour coding as in Fig. Gray 

region: experimental range of B(K + — > tt + vv) (see Eq. W. 6]) ). Black horizontal line: 
upper bound of B(Kl — > H + (see Eq. $6.18\) ). 



• Without the upper bound on BiK^ — > fi + fi ) the branching ratio B(K + — > 7T + uu) 
could be much larger. 

• The fact that the increase of one of the two branching ratios implies uniquely 
the increase of the other one signals the dominance of left-handed currents in NP 
contributions. 

Concerning the latter point, if right-handed couplings were dominating as found in RSc 
scenario [37] and some Z' scenarios in [3], we would find an anti-correlation i.e. an 
enhancement of B(Kl relative to the SM would imply suppression of B(K + — > 

Tt + vv) and vice versa. This different behaviour originates in the fact that the K + — > 7i + uu 
transition is sensitive to the vector component of the flavour violating Z coupling, while 
the Kl — > fi + n~ decay measures its axial component. In other words, the correlation 
between K + — > 7i + uu and — > offers a clear test of the handedness of NP 

flavour violating interactions and Fig.[9]shows transparently that in TUM the left-handed 
couplings are at work. A more general discussion of these points can be found in [92] • 




8.6 K L 7T°£+£- 

In the left panel of Fig. [10] we show B(Kl — > 7r°e + e~) vs B(Kl — > 7r°/i + yU— ). In the 



right panel the correlation between B(Kl — > 7i°e + e ) and B[Kl — > n^vv) is shown. 
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Figure 10: B(K L -)■ 7r°e + e~) versns £>(-K"l ->■ n°fi + fi-) (left) andB(K L -»■ 7r°z/z/) (right). 
Colour coding as in Fig. 




B(K%i*vy) [10 -11 ] 



Figure 11: -B Si d — >■ versus B(K + — )■ 7r + z/z/) and B(K L — > ifivv). Gray region: exp. 

range. Colour coding as in Fig. 



We observe that the correlations between all three branching ratios are very strong and 
all branching ratios are enhanced relative to the SM values but NP effects in B(Kl — > 
7r°£ + £~) are as expected much smaller than in B(Kl — > 7r°uu). We recall that we only 
show the results for the constructive interference between SM and NP contributions. 



8.7 -B s ,d - > versus K + — > tt + w and Kj — y Tt Q w 

In view of the small number of parameters in the TUM these decays are correlated with 
each other. In the left panel of Fig. [TT] we show the correlation between B(B S — > 
and B(K + — > -k + vv). This correlation is similar to the one between B(K^ — > -n^vv) and 
B(K + — > 7i + uu) as B(B S — > / u + ^r) similarly to BiK^ — > n°vv) depends primarily on 
. Not surprisingly the correlation between B(B S — > and B(Kl — >■ 7r°z/z/) is very 

simple. We show it in the right panel of Fig. [TTJ 
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Figure 12: BiB — > K*vv) versus B(B S — > fi + {i~). Gray region: exp. 1 and 2a range of 
B(B S — > n + (see Eq. ^6. 30\) ). Colour coding as in Fig. O 



8.8 B — ^ {X s , K, K*}vv 

As the right-handed currents are suppressed in TUM and rj f» 0, the three branching 
ratios in question are strongly correlated with each other and knowing one of them gives 
the information about the other two. Moreover, as left-handed currents dominate, there 
is a strong correlation between these decays and B s ^ — > [i + fi~. In Fig. [12] we demonstrate 
this by showing B{B — > K*uu) versus B(B S — > fi + fi~). 

8.9 Implications of b — > s£ + £~ Constraints 

Presently the NP effects found by us are consistent with the experimental data on B s ^ — > 
although with improved upper bound on B s ^ —> A i+ A i ~, the TUM could have 
problems with describing the data. However, also the data on B — >■ X s £ + £~, B — > 
K*£ + £~ and B — > Kl + l~ improved recently by much and it is of interest to see whether 
they have an impact on our results. A very extensive model independent analysis of 
the impact of the recent LHCb data on the Wilson coefficients and CfJ has been 
performed in [93] and we can use these results in our case. Other recent analyses of 
b — > s£ + £~ can be found in [9^195]. 

In view of suppressed vectorial couplings of Z to muons the bounds on Cg are easily 
satisfied. Therefore we will only check whether for the ranges of parameters considered 
by us the resulting coefficients C{g satisfy the model independent bounds in [93]. As 
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these coefficients are scale independent we can use the formulae in (I5.18P and ( 15 .20 j) and 
compare the resulting coefficients with those in the latter paper. 

We find first that in TUM 

sm 2 9 w Cl P = -AX L (B S ) = -0.31 ( ^XYef, (8.10) 



\ M J 



sin 2 e w C' w = -X R (B s ) = -0.31 -L 



2 ' r 2.0-5.H 2.2 



+ 



Y2 1 y2 

^13 ^23J 



10~ 3 . (8.11) 



On the other hand the allowed 2a ranges of Cj'J are shown in Figs. 1 and 2 of [93]. They 
are given approximately as follows: 

- 2 < Re(C[ ) < , -2.5 < Im(C%) < 2.5 , (8.12a) 

- 0.8 < Re(C 1 N p ) < 1.8 , -3 < Im(Cio) < 3 . (8.12b) 

Especially, the new data onB-f K*fj, + fi~ allow only for negative values of the real part 
oiC' w 

Re(C[ ) < 0. (8.13) 

As seen in (18. lip in TUM this condition is satisfied as Re(C[ ) is predicted to be negative. 
Moreover, in view of the suppression of right-handed currents in TUM, these data have 
no impact on our results. 

On the other hand the lower bound on the real part of Cf p gives an upper bound on 
the ratio e^/M: 

3 ~^) e? < 0.78. (8.14) 

Therefore for M = 3 TeV these data prefer close to its lower bound of 0.8. Conversely 
this bound implies 

M > 3 TeV. (8.15) 

If this behaviour will be confirmed in the future by more accurate data and improved 
form factors that enter the analysis of |93j, the TUM will favour this low value implying 
rather precise values for various branching ratios for M = 3 TeV. In particular we 
predict then: 

B(B S ->-/i + /0 = (4.6 ±0.4) x 10" 9 , B(B d ^fi + fi-) = (1.5 ±0.1) x 10~ 10 (8.16) 
and 

B(K L 7i°uu) = (3.8 ± 0.4) x 10 -11 . (8.17) 
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As these data do not fix A23 the branching ratio B(K + — > ti + vv) is still consistent with 
all the data for 

9 x KT 11 < B{K + -»■ ti + vv) < 16 x 10~ u . (8.18) 

We should remark that all the branching ratios depend to an excellent approximation 
on the ratio e^jM. Therefore the predictions given above can be kept, while satisfying 
([HUD, for 3 TeV < M < 3.8 TeV by increasing e£. As < 1, for M > 3.8 TeV the 
values of the branching ratios will decrease with increasing M. 

9 Summary and Outlook 

In the present paper we have performed a detailed analysis of particle-antiparticle mixing 
and of the most interesting rare decays of K and B mesons in the MTFM concentrating 
on its simplest version, the trivially unitary model (TUM), in which the Yukawa matrix 
in the heavy down quark sector is unitary and the corresponding matrix in the heavy up 
quark sector is a unit matrix. The new contributions to FCNC processes are dominated 
by tree-level flavour violating Z couplings to quarks. The modifications of the P^-boson 
couplings and the generated flavour violating Higgs couplings have negligible impact 
on observables considered by us. Our analysis includes complete renormalization group 
QCD effects in NP contributions at NLO in the case of AF = 2 and in rare K and B 
decays. 

The TUM can correctly reproduce the masses of quarks and the CKM matrix leaving 
for fixed heavy quark mass M three real and positive definite parameters which together 
with CKM couplings and quark masses govern NP contributions to FCNC processes. 
The paucity of free parameters in the TUM implies a number of correlations between 
various observables within the K system, within the B s d system and between K and 
B S) d systems, and also between AF = 2 and AF = 1 observables. These correlations 
allow for a clear distinction between this model and other NP scenarios. 

The main messages of our paper are as follows: 

• The simplest version of the MTFM, the TUM, is capable of describing the known 
quark mass spectrum and the elements of the CKM matrix favouring \ V u b\ ~ 0.0037 
and 7 ~ 68°. The masses of vectorial fermions are bounded to be larger than 
M ~ 3.0 TeV implying that these fermions are still in the reach of the LHC. 
Precise lower bound would require reduction of theoretical uncertainties in FCNC 
processes. 
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• NP effects in B° sd — B St d observables are very small so that AM s> d, and S^k s 
are basically identical to SM predictions. The optimal values for \V u f,\ and 7 imply 
very good agreement of AM Si d and with the present data, while S^k s ~ 0.72 
is by 2a higher than its present experimental central value. It will be interesting 
to see how future LHCb results compare with this prediction. 

• NP effects in Ek could in principle be sizable but the interplay of NP contributions 
from LL and LR operators to the data for the later and the upper bound 
on B(Kl — > /i + yU~)sD makes also NP effects in ek very small. Simultaneously 
right-handed flavour violating Z couplings to quarks are forced to be suppressed 
leaving the corresponding left-handed couplings as the dominant source of NP 
contributions to rare K and B s ^ decays. 

• The pattern of deviations from SM predictions in rare B decays is CMFV-like with 
an important prediction not common to all CMFV models: B(B s d — > f^ + ^ ) are 
uniquely enhanced by at least 35% relative to SM values and can be by almost 
a factor of two larger than in the SM. While still consistent with LHCb results, 
B(B S — > may turn out to be too high to agree with the future improved 
data. Finding this branching ratio to be enhanced would be good news for the 
TUM. Also b — > svv transitions are enhanced by a similar amount. 

• The model predicts uniquely the enhancements of the branching ratios for K + — > 
n^vv and Kl — > tt°uu by similar amount as the rare B decay branching ratios. 
In particular the correlation between Kl — > 7r°uu and B s d — > n + fi~ is CMFV-like 
but the correlation between Kl — > -n^vv and K + — > 7i + uu shows a non-CMFV 
behaviour. NP effects in Kl —> tc°£ + £~ are found to be significantly smaller. 

• The implications of the recent data on other b — > s£ + £~ transitions from LHCb is 
to suppress some of the effects listed above so that eventually TUM makes rather 
sharp predictions for all FCNC observables that in the case of rare decays differ 
from the SM: see (IQty - tETTgj) . 

In view of these very definite predictions we are looking forward to improved experimen- 
tal data and improved lattice calculations. The correlations identified in the TUM will 
allow to monitor how this simple NP scenario discussed by us face the future precision 
flavour data. In case of difficulties a possible way out would be to increase the value of 
M or take X u ^ t which in turn would allow to introduce two CP phases in X D and 
use the parametrization in (12 . 1 3[) for it. This would have an impact on CP- violating ob- 
servables with smaller effects on CP-conserving ones, although it could allow in principle 
suppressions of various rare decay branching ratios which is not possible in the TUM. 
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It would also introduce NP effects in the D system. Finally the Yukawa couplings of 
vectorial fermions could be non-unitary matrices. But these generalizations are not yet 
required. 
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A Rotation matrices from mass diagonalization 

In Section [3] we use the following shorthand notation: 



A 12 


— A 13 A 32 


A 12 A 33 


A 22 


- \ U \ U 

— A 23 A 32 


\ U \ U 
A 22 A 33 


A 13 


- \ U \ U 

— A 13 A 22 


\ U \ U 
A 12 A 23 



A 33 = |A 



U |2 
33 1 



x 21 



A 31 



A 23 A 31 

\ U \ U 
A 22 A 31 



A 21 A 33 

\ U \ U 
A 21 A 32 



(Al) 
(A.2) 
(A.3) 
(A.4) 



A 12 


= A 13 A 32 


~ A 12 A 33 


A 22 


— A 23 A 32 


A 22 A 33 


A 13 


= A 13 A 22 


X D X D 
A 12 A 23 



X D 
A 



21 

D 
31 



A 23 A 31 



A 22 A 31 



X D X D 
A 21 A 33 



A 21 A 32 



(A.5) 
(A.6) 
(A.7) 



\D — \D \D* , C D2\D \D* \D _ \D \D* , C D2\D \D* ( a q\ 

A 13 = A 13 A 33 ' £ 23 A 12 A 32 A 23 = A 23 A 33 "+" £ 23 A 22 A 32 

A 33 = l A 33| 2 + £ 23 2 | A 32| 2 (A-9) 

Here we give explicit expressions for the entries of the rotation matrices V^'j^ in Eq. (13. 9ft - 
(EIl: 

\U \U \U \U* 

A 22 A 33 A 33 A 22 

\U* \U* \U* \U 

n,R — _2A_ R _ *31. R _ ^32_ R _ *31 / A ii\ 

U l ~ ru* ' 2 ~ \U* ' M 3 - \U* ' 4 \[7 ' l A - Li J 

A 22 A 33 A 33 A 22 
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d\ 

df 



A 12 
A 22 



4 



A 13 

X D 

A 33 



21 , D2 A 32 A 31 



\D* 
A 22 



; 23 \D t D 
A 33 A 22 



4 



I X D I 
l A 33l 



A 23 

X D 

A 33 



4 



A 13 
A 22 



di 



\D* \\D\ 
A 31 \ A 33\ 
\D* 
A 33 



(A.12) 
(A.13) 



df 



ID ' °4 
A 22 



I X s I 
l A 33l 



\D* \\D\ 
A 32 \ A 33\ 
\D* 
A 33 



d R --^ 

A 22 



To make a comparison with the RS scenario in [6] easier we give here the 
between the different notation used for the rotation matrices (up to phases) 



12 



d[ 



IjJ 



IjJ. 



21 



-d{ 



13 



d L 

a 2 ) 



LO. 



and similarly with d <-> tt, and 



Pl2 
P21 



P?2 
P22 



A' 
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pr 3 



31 



2 ) 



d L 



23 



LO 



d L 3 , 



P23 



U 



32 



R 
3 • 



-d L * 
a 3 5 



— M 



ft* 



P31 



U 



R 
4 ) 



P32 



-U 



R* 



R 



d{ 



P d 33 



d 



Pl3 ~ ^2 j 

-R d jR 

4 > P31 ~~ u 6 



P21 



d 



R 
3 ) 



P32 



P23 
-df . 



d 



R 

5 ) 



(A.14) 
translation 

(A.15) 



(A.16) 



(A.17) 



B Flavour violating Z and W couplings 

Due to the mixing of SU(2) doublets and singlets flavour violating Z and Higgs cou- 
plings are induced and the W-couplings are modified. These effects are parametrized by 
hermitian matrices A^' 1 ^ and A^ D defined already in p] and readdressed in Section [3j 
Here we give analytic formulae for these couplings in the mass eigenstate basis (defined 
in Section [3]): 



1 



Ml 



(- 



rr-jT 



cc. 
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e ujj U\ 



cc. 



„icu m u mt R 
6 UJJ u 2 
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(l + \uf\ 2 ) e 



i(c v -b v ) iricint 



™L (1 + 
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R * U f) 
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(B.18) 
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CC. 
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(B.19) 
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e .Q.O u 2 

fc l fc 3 
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where "cc." denotes the complex conjugate of the related entry. 



(B.21) 
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In the unitary model without the constraint in X u we get: 
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where of course the uf ,R , df' R are not all independent because they are functions of the 
unitary \ U)D . Inserting the angle-parametrization from Sec 12.21 we get for Aj? R the same 
as in TUM but very lengthy expressions for A V LR which we do not list here. 
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C Approximate Expressions for Z contributions to 
AF = 2 Amplitudes in TUM 



The tree-level Z contributions to mixing amplitudes M^ 2 and M\ 2 are given approxi- 
mately as follows 
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In the K sector the above expressions directly give the new contribution 5sk and 5AM K , 
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In the B-sector one can obtain the relevant quantities 
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D Approximate Expressions for Z contributions to 
AF = 1 Amplitudes in TUM 

obtain approximate expressions for the relevant 



Proceeding in the same manner we 
quantities in AF = 1 observables 
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Moreover one has 

AY d = AY S = AX L (B d ) = AX L (B S ) Y' d = X R (B d ) 
AY A (K)=AX L (K)-X R (K) 



Y' s = X R (B S ) 
AZ' q = AY; AZ q =AY q (D.45) 
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